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EXECxrrivE  summary 


Tliis  report  presents  in  detail  the  the  theory  support  of  DNA’s  Plasma  Radiation  Source  (PRS) 
RAST  program  carried  out  by  NRUs  Radiation  Hydrodynamics  Branch  (Code  6720)  in  FY  1995. 
The  major  and  continuing  goals  of  the  PRS  program  are  focused  in  two  areas:  improved  yields 
and  fidelity.  Put  another  way,  both  adequate  and  accurate  doses  must  be  available  for  meaningful 
effects  testing.  Of  course,  unlimited  resources  are  not  available  to  accomplish  this  mission,  and 
one  promising  option  for  cost-effective  testing  is  embodied  in  ACE  4,  the  most  powerful  of  the 
inductive  energy  store  (lES)  machines  yet  built.  The  work  described  in  this  report’s  nine  sections 
is  aimed  not  only  at  the  general  objectives  of  yield  and  fidelity,  but  specifically,  how  to  accomplish 
them  in  the  environment  of  the  particular  challenges  posed  by  ACE  4.  The  subjects  discussed  are, 
section  by  section: 

(1)  A  series  of  initial  A1  wire  array  experiments  on  ACE  4  was  carefully  designed  using  all 
of  the  knowledge  which  has  been  obtained  by  careful  benchmark  experiments  on  other  high- 
current  machines.  ACE  4  performed  poorly  comfW  to  these  other  machines  (Double  EAGLE, 
Phoenix,  Saturn,  and  Blackjack  5),  with  K  shell  yields  about  an  order  of  magnitude  below  that 
routinely  obtained  at  similar  currents.  This  section  discusses  what  is  known  about  the  load 
dynamics  requirements  for  good  yield,  as  well  as  the  thorough  and  systematic  experiments  on 
several  machines  which  have  provided  such  knowledge.  Although  load  stability  is  undoubtedly 
an  important  issue,  the  initial  poor  results  on  ACE  4  probably  resulted  from  power  flow  problems. 
Methods  of  correcting  and/or  working  around  these  problems  are  suggested  and  will  be  pursued. 

(2)  One  of  the  requirements  for  greater  fidelity  is  to  produce  harder  x-rays.  This  means  stripping 
higher  atomic  number  elements  of  more  electrons,  which  will  inevitably  strain  available  machine 
energies.  This  section  describes  an  extensive  preliminary  investigation  of  this  problem  earned  out 
with  detailed  radiation  hydrodynamic  simulations.  Specific  results  are  presented  for  Al,  Ar,  Ti,  and 
Kr.  Two  of  the  important  conclusions  are  that  there  was  a  Z  dependent  lag  in  ionization  into  the  K 
shell  that  adversely  affects  the  emission,  becoming  worse  with  increasing  atomic  number,  and  that 
opacity  and  irmer  shell  absorption  were  also  important  limiters  of  the  K  shell  power.  The  plasma 
loads  generally  remained  in  the  inefficient  regime  in  these  simulations,  in  which  emission  scales  as 
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the  square  of  the  mass. 

(3)  One  option  for  obtaining  hard  photons  that  initiaUy  seems  very  attractive  is  to  employ  nonthermal 
processes  in  which  hot  electrons  either  produce  bremsstrahlung  within  the  load  plasma  or  excite 
inner  shell  K  alpha  radiation.  In  this  section  it  is  shown  quantitatively,  using  optimistic  assumptions 
regarding  the  internal  pinch  physics,  that  these  processes  ate  so  inefficient  that  they  are  not  likely 
to  be  interesting  sources  in  the  20-100  keV  range. 

(4)  A  possible  technique  to  improve  both  yield  and  fidelity  in  the  mass  squared  regime  of  machine¬ 
load  coupling  is  to  implode  a  gas  puff  onto  a  wire  located  on  axis.  1-D  modeling  of  the  radiative 
implosion  dynamics  of  such  an  approach  is  described  in  Sec.  IV.  Theory  indicates  that  at  small 
mass  loads  the  puff-on-wite  configuration  should  outperform  the  puff-only  loads  on  ACE  4.  The 
presence  of  the  on-axis  wire  serves  to  increase  the  ion  density.  An  initial  series  of  experiments  on 
Double  EAGLE  in  which  an  Ar  puff  stagnated  onto  a  central  Ti  wire  showed  significant  increases 
in  the  yield  above  photon  energies  of  4  keV.  This  bodes  weU  for  the  technique  as  a  fidelity  as  well 
as  yield  enhancer. 

(5)  The  Rayleigh-Taylor  (RT)  instability  is  a  potential  threat  to  the  viability  of  PRS  loads,  especially 
when  the  initial  radius,  runin  time,  and/or  current  risetimes  are  large.  The  RT  instabihty  occurs  in 
the  Z  pinch  environment  because  a  magnetic  field,  with  essentially  no  mass  or  inertia,  accelerates  a 
fluid,  namely  the  pinch  plasma,  of  finite  mass.  X-ray  yields  are  reduced  when  this  instability  grows 
too  rapidly  since  some  fraction  of  the  radial  kinetic  energy  will  be  transformed  to  axial  kinefic 
energy,  and  the  compression  will  be  disrupted.  Section  V  presents  detailed  results  for  growth  rates 
for  loads  smoothly  varied  from  a  sheU  to  a  uniform  fill.  A  tradeoff  is  noted  between  the  enhanced 
stabiUty  of  filled  cylinders  and  the  attendant  reduction  in  load  kinetic  energy.  The  calculated  yields 
of  the  stable  loads  are  superior  to  those  of  thin  shells,  indicating  in  this  case  that  the  tradeoff  is 
worthwhile.  An  investigation  of  more  complex  tailored  density  distributions  has  been  undertaken, 

and  will  be  continued  in  the  future. 

(6)  Countering  the  threat  of  RT  instabiUties  by  deUberately  inducing  rotation  in  Z  pinches  is 
investigated  in  this  section.  It  is  found  that  such  a  technique  does  help  stabilize  long-wavelength 
modes  in  the  run-in  phase.  However,  the  short-wavelength  modes  ate  not  likely  to  be  affected  by 
rotation,  and  therefore,  to  improve  the  yield  and/or  power  in  PRS  appUcations,  this  method  should 
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be  employed  in  combination  with  other  methods  of  stabilization  such  as  tailored  density  profiles. 

(7)  For  years,  the  K  shells  of  several  elements,  such  as  neon,  aluminum,  and  argon,  have  served  the 

role  of  both  providing  the  radiation  needed  for  DNA’ s  simulation  c^ability  and  diagnosing  the 
pinch.  There  is  of  course  no  fundamental  reason  why  other  atomic  shells  should  not  be  satisfactory 
for  such  purposes,  especially  when  enhanced  fidelity  is  desired  and  the  load  design  is  constrained 
by  the  chemical  or  mechanical  properties  of  the  materials.  The  L  shell  of  krypton,  for  instance, 
radiates  copiously  in  the  1.6-2.0  keV  region  of  the  x-ray  spectrum.  This  section  presents  the  results 
of  an  effort  to  bring  the  Branch’s  L  shell  data  base  closer  to  the  standard  of  comprehensiveness  set 
by  our  extensive  K  shell  capability.  Of  great  importance  is  the  discovery  of  several  new  diagnostics 
of  the  plasma  conditions,  including  the  use  of  satellite  to  resonance  line  ratios  to  infer  the  heating 

rate  of  the  pinch. 

(8)  The  ultimate  test  of  any  theory  or  numerical  simulation  is:  does  it  agree  in  large  measure 
with  experiment?  Can  it  be  used  to  interpret  experimental  results  and  also  to  study  tradeoffs 
in  experimental  design?  To  ascertain  the  degree  to  which  this  is  true  or  false  requires  that 
the  simulation  predict  quantities  which  are  measureable,  preferably  those  which  are  readily  and 
routinely  measured.  Much  of  the  data  from  Z  pinch  experiments  consists  of  filtered  spectrally  and 
spatially  resolved  radiation  profiles.  This  section  describes  progress  and  results  in  developing  the 
capability  to  predict  pinhole  images  of  pinches,  resolved  and  spectrally  filtered  to  whatever  degree 
is  requited  to  conform  to  the  capabilities  of  the  instruments  deployed  on  a  particular  experiment  An 
example  of  how  this  diagnostic  comparison  will  be  carried  out  in  practice  is  given  for  an  aluminum 

wire  array  implosion. 

(9)  Last  year’s  report  described  how  both  yield  and  fideUty  were  improved  by  using  Mg-coated  A1 
wires  on  Double  EAGLE.  Reduction  of  the  average  line  opacity  improved  the  K  shell  emissivity, 
and  aHHing  x-rays  from  a  broader  range  of  energies  enhanced  simulation  fidelity.  To  move  toward 
optimization  of  this  highly  promising  tqjproach,  the  final  section  of  this  report  investigates  all  four 
elements  with  atomic  numbers  11-14.  It  is  concluded  that  the  method  is  amenable  to  increased 
atomic  number  elements,  and  that  had  the  Mg  coat  been  a  Si  coat,  the  yield  would  have  doubled 
instead  of  increasing  by  28-50%.  Hiis  conclusion  assumes  that  the  Si  cold  start  would  be  similar 
spatially  and  thermally  to  that  of  the  Mg  in  the  Double  EAGLE  experiments. 
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I.  DESIGNING  AND  DIAGNOSING  ACE  4  ALUMINUM 
ARRAY  EXPERIMENTS 


A.  Backgrounds  Double  Eagle,  Saturn,  BJ5,  and  Phoenix  Experiments 

ACE  4  is  the  first  inductive-store  generator  to  be  considered  for  use  as  a  plasma  radiation 
source  (PRS)  for  NWET  ^plications.  Some  PRS  work  had  been  done  on  the  FALCON  generator 
at  Physics  International  Co.  (PI)  to  demonstrate  the  role  of  a  plasma  opening  switch  (POS)  on  PRS 
performance.  It  was  used  to  shorten  the  current  risetime  and  to  increase  K-shell  yields  from  given 
initial  diameter,  neon  gas  puff  implosions.^  The  shortened  risetime  allowed  more  kinetic  energy  to 
be  generated  starting  from  the  given  initial  diameter.  The  larger  implosion  energy  allowed  neon  to 
enter  a  region  in  (gas  puff  mass)  -  (implosion  velocity)  space  where  the  most  efficient  conversion  of 
kinetic  energy  into  K-shell  radiation  takes  place.  Generally,  an  inductive-store  machine  is  designed 
to  have  a  long  current  risetime  and  a  POS  that  shortens  this  risetime.  ACE  4  was  designed  as  a 
bremsstrahlung  machine  with  a  large  inductance  and  with  this  long  risetime  feature  built  in. 

Because  ACE  4  has  a  relatively  long  current  risetime,  with  or  without  an  opening  switch,  it 
is  necessary  to  do  large  radius  PRS  experiments  in  order  to  ignite  the  K-shell  of  either  aluminum 
or  argon.  In  general,  as  pulse-power  machines  scale  up  in  energy,  it  is  necessary  to  initiate  PRS 
implosions  from  larger  radii.  For  this  reason,  the  viability  of  large  radii  implosions  needs  to  be 
assessed  on  existing  non-inductive-store  machines  in  order  to  establish  benchmarks  for  similar 
experiments  on  ACE  4. 

Aluminum  Double  Eagle  experiments  were  designed  late  last  fiscal  year  to  increase  K-shell 
yields  at  large  array  radii  (and,  therefore,  for  longer  run-in  time  implosions).  These  experiments 
achieved  higher  K-shell  yields  from  larger  array  diameters  than  had  ever  been  achieved  before.^ 
The  measured  yields  and  the  calculated  conversion  efficiencies  (of  kinetic  energy  into  kilovolt  x 
rays)  are  listed  in  Table  I.  The  newly  measured  Double  Eagle  yields  differed  in  another  way  from 
the  old  yields:  they  had  a  maTimiiTn  that  was  shifted  to  higher  implosion  velocities  (see  Fig.  (1)). 
This  put  Double  Eagle  in  closer  agreement  with  SaUim  and  Phoenix  results  than  it  had  been  before.® 
These  three  sets  of  experiments  clearly  established  that  K-shell  yields  were  maximized  in  a  region 
of  (load  mass)  -  (implosion  velocity)  space,  shown  in  Fig.  (1),  that  is  predictable  fr-om  present 
theory.  This  result  provides  an  important  benchmark  for  ACE  4. 

The  analysis  and  interpretation  of  some  Saturn  A1  spectral  data  was  also  finished  late  last  year 
and  documented.'*  The  data  was  several  years  old  and  had  been  obtained  through  partial  DNA 
support.  It  consisted  of  pinhole  photographs,  which  were  scanned;  PCD  traces  giving  the  x-ray 
pulsewidth  of  the  implosion;  digitized  spectral  data;  and  electrical  and  thermal  measurements  of 
power  and  yield.  The  K-shell  yield  data,  which  is  shown  in  Fig.  (2)  of  Ref.  (5)  on  page  42,  exhibited 
a  well-defined  peak  when  plotted  as  a  function  of  the  diameter  of  the  wires  used  in  each  of  the 
loads.  Software  was  developed  to  automatically  determine  a  consistent  set  of  pinch  parameters 
from  the  complete  set  of  spectral  data.  These  parameters  included  pinch  diameter,  pulsewidth, 
average  power  output,  and  spectral  line  ratios.  The  observed  power  and  the  computed  fine  ratios 
were  then  compared  to  theoretical  predictions  of  these  quantities,  which  had  been  calculated  as 
a  function  of  an  average  electron  temperature  and  an  average  ion  density  of  the  observed-sized 
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The  date  analysis  that  produced  Fig.  (3)  of  Ref.  (5)  provided  an  essential  insight  into  te 
dynamics  of  wire  array  z-pinches.  It  suggested  thal  when  the  wues  were  too 
explode  properly,  and,  when  they  were  too  small,  they  did  not  implode 
trend  in  radiating  mass  proved  that  wire  thickness  is  an  important  factor  m  determining  K 
S  ^ially  on  a  short  risetime  machine  like  Saturn.  A  new  set  of  large  ^ay  dr^eter^ 
Lperim^  WK  then  designed  to  demonstrate  that  one  could  woit  around  the  neganve  yield 
fintogs  that  had  shown  up  in  the  early  Saturn  experiments.  new 

to  Sandia  during  the  DNAffrench  data  exchange  meetmg,  which  was  held  at  Sania.  They  were 
S^luy  performed  on  Saturn  shortly  thereafter  and  their  yield  behavior  was  discussed  at  the 


DNA  PRS  Workshop  by  C.  Deeney. 

The  shot  matrix  for  this  second  set  of  Saturn  experiments  is  given  in  Table  ^  The  location  of 
these  shots  in  (load  mass)  -  (implosion  velocity)  space,  i.e.,  in  (m  -  Vf)  space,  is  shown  in  ig^(  )• 
Sy -he  ex^rinients  -re  Ligned  to  investigate  yield  behavior  around  tiie  upper  bouridt^ 
(niBpWi,  onhe  efBcient  emission  region.  This  (mass  breakpomt)  boun^  was 
using  soft  hnplosion,  x-pinch  calculations.'  The  location  of  the  boundary.  (mBp)H.,^.  that  is 
3Ld  from  hard  implosion  calculations^  is  also  shown  in  Fig,  (2).  As  the  30  wire  dampoint 
«  26^cm  shows,  Jre  was  tittle,  if  any,  fatioff  in  yield  near  the  hard  implosion  bouuda^  m 
these  experiments.  This  demonstrates  that  some  of  the  softemng  of  implosions  may  be  caused  by 
SX^nirbulence  when  an  insufficient  number  of  wires  are  used  m  the  implosion. 

These  experiments  were  designed  to  increase  the  mass  participation  in  K-shell  emission  on- 
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axis  by  increasing  the  number  of  wires  in  the  arrays,  by  decreasing  their  size,  and  by 
the  return  current  path  out  to  a  radius  of  2.5  cm  in  order  to  reduce  implosion  asymmetnes.  The 
larger  wire  number  and  larger  return  current  radius  increased  the  symmetry  and,  therefore,  the  1-D 
character  of  the  implosions.  The  smaller  wire  sizes  allowed  the  wires  to  explode  more  completely 
before  implosion  began.  However,  by  increasing  the  wire  number,  the  current  flow  to  e^h  wire 
was  reduced.  The  low  yields  of  the  1.2  and  1.0  mil  42  wire  shots  indicate  that  this  reduction  may 
have  contributed  to  an  inadequate  explosion  of  these  wires.  When  only  30  wires  were  used,  the 
1 0  mil  wire  shot  performance  improved  markedly.  The  30,  1  mil,  wire  shot  also  outperformed 
the  best  yielding  of  the  shots  in  Fig.  (2)  of  Ref.  (5).  This  shift  up  the  mass  breakpornt  curve  of 
the  tnayimnm  yielding  of  the  aluminum  wire  shots  on  Saturn,  which  was  predicted,  is  displayed  rn 

Fig.(l). 

These  experiments  also  produced  a  somewhat  unexpected  result.  While  the  yields  of  the 
42  wire  shots  were  sUghtly  lower  than  those  of  the  nearby  30  wire  shots,  the  x-ray  pulsewidths 
of  the  42  wire  shots  were  shorter  than  those  of  the  30  wire  shots,  i.e.,  the  x-ray  powers  had 
increased.  This  suggested  that  both  x-ray  yield  and  pulsewidth  could  be  independently  controlled 
by  carefully  adjusting  both  the  wire  number  and  wire  size  in  PRS  loads.  To  test  this  idea,  a  new 
set  of  experiments  needed  to  be  done  in  which  the  array  mass  was  fixed  and  the  number  of  wires 
was  varied.  The  design  of  such  experiments  is  facilitated  by  plotting  mass  contours  in  (wire-size)  - 
(wire-number)  space.  Two  were  computed.  They  are  shown  in  Figs.  (3)  and  (4)  for  A1  and  Ti  wires 
respectively.  The  circles  on  the  A1  wire  contours  show,  for  example,  how  one  could  approximately 
construct  a  200  load  from  0.6, 0.7, 0.8, 0.9,  or  1.0  mil  wires  by  using  40,  30, 23, 18,  or  14 
wires  respectively.  These  contours  were  presented  to  Sandia,  and  they  stimulated  a  new  series  of 
experiments  on  Saturn  that  did,  indeed,  produce  x-ray  pulsewidth  variations  for  approximately  a 
given  K-shell  yield.  However,  the  total  yields  did  go  up  as  the  power  increased.  As  the  wire  number 
increased  and  the  gap  spacing  between  wires  decreased,  the  implosions  became  more  symmetric 
and  1-D-like.  These  experiments  indicate  that  arrays  will  implode  to  the  same  final  pinch  radius, 
independent  of  the  initial  array  radius,  when  the  gap  spacing  between  wires  is  initially  less  than  1 
mm.  These  findings  will  need  to  be  carried  over  and  tested  on  ACE  4. 

Finally,  because  Double  Eagle,  Phoenix,  and  Blackjack  5  (BJ5)  generate  comparable  currents, 
it  is  puzzling  that  the  maximum  yields  seen  on  Double  Eagle  and  Phoenix  occurred  at  different 
locations  in  the  efficient  emission  region^’®  (see  Fig.  (1)).  Moreover,  the  Double  Eagle  yields 
appear  to  be  in  agreement  with  earUer  published,  BJ5  yields.**®  However,  this  BJ5  data  w^ 
incomplete.  More  data  is  needed  to  help  resolve  the  difference  between  Phoemx  and  Double  Eagle 
results  and  between  Double  Eagle  and  Saturn  results.  However,  during  the  year,  BJ5  expenments 
had  to  be  piggybacked  onto  debris  and  diagnostic  testing  at  MU.  This  testing  placed  contraints  on 
the  kind  and  number  of  experiments  that  could  be  done,  so  the  results  obtained  from  them  were 
suggestive,  but  not  definitive.  Two  centimeter  length  wires  could  not  be  used,  for  example,  nor 
could  small  diameter  arrays  or  large  wire  number  arrays.  Undiagnosed  changes  in  machine  power 
flow  had  also  occurred,  so  that  BJ5  could  not  reproduce  earUer  yield  benchmarks  that  had  been 
obtained  with  a  canonical  MLI  load.  The  yields  were  lower  as  a  result  of  power  flow  degradaUon. 

The  proposed  BJ5  experiments  and  the  measured  yields  are  shown  in  Table  HI.  The  expenments 
were  designed  to  explore  regions  oim-vf  space  where  Phoenix  and  Double  Eagle  had  their  best 
performance  (see  Fig.  (1)),  albeit  with  2  cm  length  loads.  There  was  also  an  attempt  to  compare 
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length  arrays  had  produced  low  yields  for  the  3  cm  shots  that  could  not  be  explained,  "nieir 
x-ray  conversion  efficiencies  were  less  than  30%.  Earlier  shots  on  Phoemx  with  4  and  5  cm 
length  wires  had  30%  or  more  conversion  efficiencies.  A  series  of  experiments  was,  therefore, 
planned  on  Phoenix  to  investigate  systematically  whether  or  not  an  array  length  dependence  to 
PRS  performance  existed.  DNA/RAST  had  funded  these  experiments  a  year  ago  on  Phoenix,  but 
they  were  not  done  because  of  machine  problems.  The  power  flow  on  Phoenix  has  smce  ^n 
improved,  transmission  Une  parameters  have  been  measured,  and  a  more  reliable  circuit  model  has 
been  derived  than  had  existed  as  part  of  the  original  machine  design.  Thus,  a  new  attempt  was 
made  to  reevaluate  the  design  and  to  carry  out  the  length  variation  experiments  this  year,  but  it  also 
failed.  Thus,  no  guidance  is  yet  available  on  how  to  optimize  load  performance  on  ACE  4  as  a 
function  of  array  length. 

B.  ACE  4  Circuit  Models 

The  biggest  risk  to  scaling  PRS  loads  to  high  current,  inductive-store  machines  is  often 
assumed  to  be  the  requirement  for  stable,  large  array  diameter  implosions.  This  assumption  is 
misleading.  An  equal,  perhaps  greater,  risk  is  that  the  power  flow  to  the  load  will  not  be  well 
understood  or  controlled.  Load  performance  degradation  occurs  at  any  array  diameter  when  power 
flow  symmetry  is  not  maintained  since  asymmetries  drive  load  instablities,  which  lead  to  ba 
implosions,  and  bad  implosions  do  not  have  large  x-ray  yields.  The  Double  Eagle  and  Saturn 
experiments  described  above  demonstrate  that  pulse  power  generators  can  be  built  with  synmetnc 
and  reproducible  power  flow,  which  can  be  described  by  dependable  circuit  models.  Load  designs 
can  then  be  found  that  convert  6%  and  more  of  theMarx  bank  energy  into  kilovolt  x  rays. 

ACE  4  was  originally  conceived  to  be  a  relatively  low  (26  nH)  inductive  store  coupled  by 
way  of  a  low  inductance  MTTL  (10  nH)  to  PRS  loads  (see  Fig.  (6)).  The  short-circuit  load  current 
from  this  generator  is  shown  in  Fig.  (7).  It  rises  to  8  MA  in  roughly  120  ns  foUowing  the  opemng 
of  the  POS  at  1  fis.  This  version  of  ACE  4  is  capable  of  imploding  a  4  cm  long  aluminum  load 
of  328  fig/cm  from  a  radius  of  3  cm  to  energies  in  excess  of  300  kJ.  'Diis  implosion  is 
for  producing  K-shell  x  rays,  and  it  corresponds  to  roughly  7%  conversion  of  the  stored  4.66  MJ 
of  Marx  bank  energy  into  load  kinetic  energy.  One  would  expect  in  excess  of  150  kJ  of  kilovo  t 
radiation  would  be  produced  if  half  of  this  kinetic  energy  could  be  radiated  in  the  K-shell. 

The  Fig.  (6)  version  of  ACE  4  has  yet  to  be  realized.  A  later,  more  inductive  version,  shown 
in  Fig  (8)  began  to  be  tested  this  year.  Only  half  of  the  potentially  available  Marx  bank  enerp 
was  at  present,  utilized,  i.e.,  connected  to  the  PRS  load;  and  half  of  the  power  flow  through  the 
transmission  line  was  inductively  grounded  (600  nH).  The  plasma  opening  switch  was  changed  to 
a  plasma  flow  switch  (depicted  as  a  coaxial  section  with  a  variable  inductance  initially  at  21 
and  the  machine  was  operated  in  a  tandem  puff  mode.  Note  that  the  circuit  in  Fig.  (8)  has  150  ^ 
more  inductance  connecting  the  Marx  bank  to  the  load  than  the  circuit  in  Fig.  (6).  Some  of  this 
inductance  has  recentiy  been  removed.  The  94  nH  feed  was  reduced  to  60  nH;  also,  the  mductive 

ground  was  removed. 

How  the  plasma  flow  switch  works  is  uncertain.  Two  versions  of  its  operation  were  used  to 
nroiect  ACE  4  PRS  load  performance.  In  one,  the  switch  is  given  a  low  impedance  and  the  current 
into  a  short  circuit  rises  slowly  in  100  ns.  This  power  flow  model  is  labeled  2a.  If  the  switch  is  given 
a  much  higher  impedance,  the  current  into  a  short  circuit  load  rises  much  more  rapidly  (withm  20 
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ns)  This  power  flow  model  is  labeled  2b.  A  second  version  of  model  2b,  labeled  2c,  was  also 
investigated.  In  this  case  the  inductive  ground  was  eliminated  and  the  94  nH  inductance  betw^n 
the  Marx  bank  and  the  coaxial  switch  was  reduced  to  60  nH.  This  had  the  effect  of  mcreasmg 
current  flow  into  a  short-circuit  load  by  1  MA.  The  short-circuit  currents  for  each  of  these  three 
power  flow  models  are  presented  in  Fig.  (9).  For  each  of  these  cases,  the  imtial 
was  2  MJ,  corresponding  roughly  to  an  erected  voltage  on  the  Marx  capacitor  of  470  . 


Given  the  three  circuit  models  for  the  ACE  4  generator,  a  set  of  large  radius  experiments  could 
be  designed  to  evaluate  ACE  4’s  aluminum  wire  PRS  performance  capabiUties.  These  experiments 
would  indirectly  test  ACE  4’s  capability  to  transfer  current  effectively  from  the  plasma  flow  switch 
into  the  diode  region  and  then  to  a  wire  load.  A  number  of  questions  were  to  be  addressed  by  these 

experiments: 

®  Assuming  some  current  is  switched  from  the  flashboard  plasma  to  an  aluminum  array  load, 
will  the  array  implode? 

9  Is  the  tandem  puff  opening  switch  operation  understood,  i.e.,  are  the  circuit  models  correct? 
Should  the  switch  be  removed  altogether? 

®  Will  aluminum  implode  symmetrically  and  efficiently  to  produce  K-shell  x  rays  from  the  l^ge 
array  diameters  required  by  ACE  4’s  present  cfesign?  If  not,  how  much  of  a  deletenous  effect 
will  there  be  and  can  we  work  around  it? 

The  shot  matrix  for  ACE  4,  A1  wire,  PRS  experiments  that  was  proposed,  which  is  given  in 
the  following  Tables  IV- VI,  was  designed  to  address  these  questions  and  to  reveal  the  ability  oft  e 
present  ACE  4  configuration  to  transfer  current  to  the  wires  and  to  implode  them  serm-umfo  y. 
It  was  proposed  that  a  set  of  experiments  be  run  for  fixed  wire  configurations  by  v^mg  the  (tang 
Marx  voltage  for  each  configuration  between  480  and  640  kV.  The  circuit  modeling  of  ACE  4 
showed  that  it  was  far  from  an  ideal  generator  in  its  present  configuration  since,  to  produce  K-sheli 
radiation,  the  experiments  would  have  to  be  stretching  the  envelope  of  what  has  beeii  achieved  to 
date  at  large  diameters  with  wire  loads  on  other  pulse-power  generators.  The  problem  was  that 
ACE  4’s  current  risetime  required  large  diameter  loads  to  make  PRS  work,  but  that  it  did  not  have 
the  current  to  implode  enough  mass  from  large  radii.  For  this  reason,  it  was  suggested  that  4  cm 
diameter  loads  comprised  of  0.8  and  0.9  mil  diameter  wires  be  used  imtially  to  do  the  voltage  scam 
Then  if  these  shots  proved  to  be  reasonably  successful,  the  larger  6  cm  diameter  airays  of  0.8  and 
0.9  ntil  wires  could  be  used  for  a  new  scan.  An  array  diameter  of  7.4  cm  was  needed  to  do  a  useful 

1.0  mil  wire  scan. 

Six  different  experiments  in  a  family  of  experiments  using  0.6, 0.7, 0.8,  and  0.9  mil  diameter 
wires  were  carried  out.  The  arrays  were  each  of  length  £  =  3  cm  long.  They  were  compos^  of 
24  wires  with  an  array  diameter  of  4  cm.  The  yield  data  are  best  displayed  in  an  m  v  f  plot  of 
the  experiments,  which  is  presented  in  Fig.  (10).  The  2  sets  of  (circle  and  triangle)  points  l^ates 
each  ofthese  arrays  in  the  figure  by  their  initial  mass  per  length.  The  taangle  points  are  obtained 
when  the  2a  circuit  model,  described  above,  is  used  to  inwardly  accelerate  the  loads,  while  the 


l^ble  rV:  Proposed  ACE  4  Aluminum  Experiments 


Return  current  radius  =  7.0  cm,  array  length  =  3.0_cm 

Wire  Diam 

(mil) 

Wire  Mass 

(jtg/cm) 

Array  Diam 

(mm) 

#  Wires 

V 

Marx  Voltage 

(kV) 

0.6 

118 

40 

24 

1.9 

585 

0.8 

210 

40 

24 

0.95 

470 

0.8 

210 

40 

24 

1.35 

600 

0.9 

266 

40 

24 

0.95 

525 

0.9 

266 

40 

24 

1.15 

595 

Return  curre 

Ibble  V;  Proposed  ACE  4  Aluminum  Experiments 
nt  radius  =  7.0  cm,  array  length  =  3.0  cm 

_ _ _ _ _ — . — - - - - - 

Wire  Diam 
(mil) 

Wire  Mass 

(ug/cm) 

Marx  Voltage 
(kV) 

0.8 

60  ^ 

24 

490 

0.8 

24 

1.65 

635 

0.9 

266 

24 

625 

Return  curr< 

ftible  VI;  Proposed  ACE  4  Aluminum  Experiments 
mt  radius  =  7.0  cm,  array  length  =  3.0  cm 

Wire  Diam 

(mil) 

Wire  Mass 

(pg/cm) 

Array  Diam 

(mm) 

^  Wires 

V 

Marx  Voltage 
(kV) 

0.9 

266 

74 

24 

1.45 

595 

1.0 

328  ^ 

74 

24 

1.18 

585 

1.0 

328 

74 

24 

1.45 

665 

circle  points  correspond  to  the  2c  circuit.  A  different  final  implosion  velocity,  v/,  is  calculated 
depending  on  which  circuit  model  is  used.  The  maximum  measured  K-shell  yields,  taken  from 
Table  Vn,  are  listed  above  the  circle  points.  The  predicted  total  kinetic  energy  generated  in  each 
(m,  vf)  experiment  can  be  computed  for  each  case.  It  is  l/2mti^  and  is  also  listed  in  Table  VE. 
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ftimtsi 


Return  current  radius  =7.0  cm,  array  ten 


drawn  from  the  above  circmi  moaei  prcuicuuus  tuiui  - - -  .  •  i  • 

or  not  the  circuit  model  is  correct,  accurate  measurements  of  the  wire  explosion  and  implosm 
times  must  be  made.  For  this  task,  fiduciated  streak  camera  and  x-ray  diode  measurements  must 

made.  The  x-ray  yield  data  suggests  that  the  implosion  was  poor,  but  how  “ 
modei  predictions  remains  to  be  determined.  The  abiiity  of  the  current  to  be  brought  down  ftOT  7 
cm  to  2^cm  before  it  couples  to  the  wires  must  also  be  tested.  If  unmodelable  c™t  asymmetries 
are  present  in  tandem  puff  operation,  they  must  be  measured  and  then  eliminated.  There  is  already 

cxpcriiB6iit3l  evidence  for  this  problcnio 

A  number  of  different  loads  designs  can  be  tried  in  order  to  work  around  some  of  the  power 

llowpriems.  For  example,  because  of  the  -  W  - /*)  ^g  of  K-shell  emissmn  m 

die  ineffleient  scaling  t^ime,  argon  load  desi^s  that  promote 

potential  of  significantly  lowering  energy  requirements  for  obtaiiung  higher  K-sheU  yields.  The 
puff  on  wire  configuration  is  a  design  that  h; 
from  gas-puff  plasmas  because  it  promote 

nirrmT%  rrQC  nff  onto  a  tl 


regime  the  presence  of  a  relatively  low  mass  wire  on  axis  does  produce  substantially  higher  total 
K-shell  yields  than  just  the  argon  puff  gas  configuration.  Calculations  show  that  a  100  /ig/cm  argon 
gas  puff  stagnating  upon  a  single  0.8  or  1.0  mil  titanium  wire  should  give  in  excess  of  1  kJ/cm 
of  additional  K-shell  emission.  The  source  of  this  additional  emission  is  the  titamum  wire.  If, 
instead  of  pure  Ti,  a  coating  of  KCl  were  put  onto  a  carbon  or  titanium  fiber  and  then  the  imploded 
argon  would  contact  high  density  potassium  chloride  on  axis,  and  one  could  expect  both  yield  and 
spectrum  enhancements  in  the  vicinity  of  the  argon  K-series  lines. 

In  addition  to  increasing  kilovolt  x-ray  yields,  the  use  of  multi-material  load  designs  provides 
opportunities  for  diagnosing  load  performance  and  for  identifying  and  quantifying  power  flow 
problems.  One  way  to  evaluate  the  quality  of  an  implosion  is  to  gauge  its  convergence  on  axis, 
which  could  be  measured  by  its  interaction  with  an  on-axis  wire.  For  this  purpose,  spectral 
diagnostics  are  needed  that  can  distinquish  core  plasma  from  surface  plasma  conditions  in  a  pinch. 
A  small  coating  of  selenium  to  the  on-axis  wire  or  a  small  admixture  of  krypton  to  the  argon 
gas  could  provide  such  a  diagnostic.  This  L-shell  diagnostic  development  is  discussed  in  another 

section  of  this  report. 
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Figure  1.  The  experiments  on  Phoenix.  Double  Eagle,  and  Saturn  generators,  labeled  P,  D,  and 
S  respectively,  that  produced  the  most  K-sheU  emission  from  aluminum  implosions  are  shown. 
They  are  located  by  the  initial  mass  and  by  the  calculated  final  implosion  velocity,  v/,  of  the  array 
in  each  case.  The  arrows  connect  the  locations  of  the  previous  maximum  yield  experiments  on 
Double  Eagle  and  Saturn  to  the  current  ones.  The  region  predicted  to  efficienUy  convert  array 
kinetic  energy  into  K-shell  x  rays  lies  between  the  i]=l.5  and  the  {mBp)soft  boundary  curves. 
Early  predicted  yield  contours  for  3, 10,  and  30  kJ/cm  are  also  shown. 


Figure  2.  Recent  Saturn  experiments  that  probed  the  yield  behavior  of  30  and  42  wire  aluminum 
arrays  near  and  away  from  the  soft-implosion  mass  breakpoint  boundary,  {mBp)soft,  are  located 
by  their  mass  and  calculated  final  implosion  velocity.  The  measured  K-shell  yields  are  displayed 
next  to  the  experimental  locations.  Note,  there  is  a  tendency  for  the  yields  to  decrease  away  from 
the  {mBp)soft  boundary,  but  the  falloff  as  the  {mBp)Hard  boundary  is  approached  is  small  for 
the  30  wire  case  shown. 
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Wire  number 


Array  Mass  Contours 


0.6  0.7  0.8  0.9  1.0  1.1  1.2 

wire  si2e(mil) 


Figure  3.  Contours  for  given  mass/centimeter  aluminum  arrays  are  shown  as  a  function  of  the 
size  of  the  wires  in  mils  and  the  number  of  wires  used  in  the  arrays.  The  dots  on  the  200  /ig/cm 
contour  correspond  to  the  use  of  0.6,  0.7,  0.8,  0.9,  and  1.0  mil  wires  and  require  roughly  40,  30, 
23,  18,  and  14  wires  respectively. 
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Wire  size(mil) 

Figure  4.  Contours  for  given  mass/centimeter  (in  units  of  /ig/cm)  titanium  arrays  are  shown  as  a 
function  of  the  size  of  the  wires  in  mils  and  the  number  of  wires  used  in  the  arrays. 
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Figure  5.  The  locations  and  measured  K-shell  yields  of  four  Blackjack  5  experiments  are  shown. 
They  suggest  that  Blackjack  5  would  confirm  the  Double  Eagle  more  than  the  Phoenix  or  Saturn 


ACE  4  Circuit  Model  #1 


Marx  inductive 

bank  store  POS  MITL 


OAto  1.3  A 
over  80  nsec 
at  1  psec. 


Figure  6.  This  is  the  circuit  of  the  ACE  4  generator  as  it  was  originally  conceived. 
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in  roughly  100  ns  once  the  POS  opened  as  shown  in  this  figure.  This  calculated  result  assumed 
MTTL  inductance  of  15  nH. 
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ACE  4  Circuit  Model  #2 


POS  #2:  Tandem  Puff 

Figure  8.  This  is  the  circuit  of  the  ACE  4  generator  as  it  is  currently  conceived  to  operate. 
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Figure  9.  Three  calculated  currents  into  a  short-circuit  load  are  shown  for  the  tandem  puff  circuit 
model.  Curves  2a  and  2b  represent  the  current  under  the  assumption  of  a  low  or  high  impedance 
plasma  flow  switch  respectively.  Curve  2c  shows  what  happens  in  the  2b  case  when  the  inductive 
ground  is  eliminated  and  the  inductance  between  the  Marx  bank  and  the  coaxial  switch  is  reduced 

to  60  nH. 


Figure  10.  The  ACE  4  experiments  of  Table  VII  are  located  in  (m,v/)  space  differently 
depending  on  whether  circuit  2a  (the  triangles)  or  circuit  2c  (the  circles)  is  used.  The  milages  of  the 
wires  used  in  the  experiments  are  listed  under  the  triangle  points;  the  maximum  measured  K-shell 
yields  are  listed  above  the  circle  points.  Based  on  the  behavior  of  the  other  generators  (Phoenix, 
Double  Eagle,  and  Saturn),  one  might  infer  that  circuit  2a  better  describes  the  array  implosions 
than  circuit  2c.  However,  this  interpretation  is  highly  uncertain  since  power  flow  problems  were 
detected  after  the  experiments  were  finished.  It  is  also  suspected,  based  on  the  low  yield  measured, 
that  the  0.7  mil  wires  were  not  aluminum. 
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II.  INVESTIGATION  OF  K-SHELL  EMISSION  FROM 
MODERATE  ATOMIC  NUMBER  IMPLOSIONS 


A.  Introduction  .  .  .  f 

A  previous  theoretical  phenomenological  investigation  was  earned  out  mto  the  ettects  ot 

magnetohydrodynamic  (MHD)  turbulence  on  the  stagnation  dynamics  of  alummum  wire  array  and 
argon  gas  puff  z-pinch  implosions.^-®  It  showed  that  ad  hoc  increases  m  the  plasma  viscosity, 
heat  conductivity,  and  electrical  conductivity  (as  modeled  by  using  a  fixed  set  of  niultipUers  for 
these  quantities  in  one  dimensional  (1-D)  MHD  calculations)  reasonably  duplicated  the  average 
electron  temperatures,  ion  densities,  and  mass  of  the  K-sheU  emission  region  that  were  measured 
at  stagnation  for  a  variety  of  Physics  International  (PI)  aluminum  wire  and  argon  gas  puff 
experiments  While  this  phenomenological  modeling  does  not  represent  a  first  prmciples  approach 
to  incorporating  turbulence  effects  in  1-D  MHD  calculations,  it  does  suggest  the  ei^tence  of  these 
strong  multidimensional  effects,  some  of  which  may  be  produced  by  the  unequal  forces  actog  on 
small  wire  number  array  loads.  In  particular,  the  major  effect  of  increasing  the  transport  coefficients 
was  to  “soften”  the  1-D  implosions  by  decreasing  the  ion  densities  that  were  achieved  at  stagna  on 


by  several  orders  of  magnitude. 

This  “soft”  implosion  modeling  was  used  in  Ref.  1  to  reexamine  the  scaling  of  K-shell  yield 
with  load  mass  for  a  fixed  implosion  velocity.  The  results  of  this  study  suggested  that  it  would 
require  a  factor  of  6  more  mass  and  kinetic  energy  than  was  calculated  ongmaUy  usmg  h^d 
(no  enhanced  transport)  1-D  implosion  calculations;  to  attain  efficient  K-shell  X-ray  yields  (K-sheU 
yield  >  30%  of  the  kinetic  energy  that  is  coupled  to  the  load  before  stagnation)  from  a  moderate-z, 
plasma.  The  original  K-sheU  yield  scaling  predictions  were  then  modified  to  reflect  this  shift  in 
mass  and  they  were  then  compared  in  Ref.  4  with  aluminum  results  obtained  from  expe^ente 
performed  on  Saturn,  Double  EAGLE,  and  Phoenix  pulse  power  machines.  In  general  K-shell 
yield  comparisons  between  measured  Saturn  and  Phoenix  yields  and  theoretical  predictions  were, 
favorable,  whereas  Double  EAGLE  produced  over  twice  the  predicted  K-sheU  yield.  Sii^arly, 
Double  EAGLE  outperformed  both  Saturn  and  Phoenix  in  terms  of  conversion  efficiency,  which  is 
the  K-shell  yield  produced  per  unit  of  calculated  kinetic  energy  delivered  to  the  load.  The  fact  that 
Double  EAGLE  produced  approximately  twice  the  predicted  yield  points  out  one  of  the  hmitations 
of  the  scaling  calculations  to  date,  which  is  that  the  current  in  the  calculations  is  tenmnated  just 
before  stagnation  in  order  to  focus  on  the  K-sheU  emission  that  is  produced  pnmanly  due  to 
thermaHzation  of  kinetic  energy.  Since  the  load  current  does  not  shut  off  m  the  experiments,  energy 
continues  to  flow  into  the  plasma  and  to  confine  the  pinch.  Both  of  these  mechanisms  could  lead 
to  an  increased  radiative  output  that  is  not  reflected  in  previous  K-shell  jneld  seating  studies.  Mter 
experiments  suggest  why  aluminum  conversion  efficiencies  have  been  higher  on  Double  EAGLE 
than  on  Phoenix  and  Saturn.®  They  provide  some  evidence  that  the  larger  yields  may  plated  to 
a  greater  late-time  current  energy  coupling  and  plasma  confinement  achieved  on  Double  EAGLE. 


Soft  implosion  calculations  show,  because  of  the  large  kinetic  energy  needed  for  efficiently 
generating  K-shell  photons,  that  higher  Z  (Z>22)  experiments  performed  on  existing  pulse  power 
machines  wiU  optimize  K-shell  emission  at  low-f?  values,  where  i]  is  defined  as  the  ratio  of  the 
maximum  kinetic  energy  per  ion  achievable  prior  to  stagnation  to  the  minimum  energy  per  ion, 

E  ■  neededtoinstantlyheatandtoionizeaplasmaintotheKshelluponstagnation.  Forelements 

up”to  krypton  E  ■  =l  012X®-®®^  eV/ion.®  Note,  this  formula  accounts  only  for  ionization  and 
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thermal  energies  and  not  radiation  or  turbulence  losses  encountered  as  the  plasma  ionizes  to  the 
K  shell.  Since  there  have  been  no  systematic  theoretical  or  experiment  investigations  of  the 
Z  scaling  of  K-shell  emission  generated  from  low-i/  implosions,  in  anticipation  of  such 
experiments  we  report  on  the  results  of  our  theoretical  investigation  into  this  regime.  In  particular 
aluminum,  argon,  titanium,  and  krypton  plasmas  are  studied  using  a  soft  implosion  model  analysis. 
The  results  of  this  study  have  identified  the  foUowing  important  effects:  low-17  plasmas  apparently 
remain  in  an  inefficient  mass^  (or  I*,  where  I  is  current)  scaling  regime  at  much  higher  masses 
than  predicted  by  the  modified  hard  implosion  model  scaling  formulas,*  innershell  absorption  of 
photons  emitted  from  the  hot  plasma  core  and  absorbed  in  colder  outer  regions  of  the  plasma  can 
substantially  reduce  the  K-sheU  yield  for  large  mass  loads,  and  there  is  a  time  dependent  lag  in 
ionization  into  the  K-sheU  that  increases  with  the  Z  of  the  plasma.  TTiese  issues  and  how  they 
influence  the  design  of  future  experiments  are  discussed  in  the  results  and  conclusion  sections  of 

this  paper. 

B.  Model 

The  1-D  MHD  calculations  are  similar  to  those  described  in  references  6,  l,^d  7.  They 
contain  a  comprehensive  time  dependent  description  of  the  ionization  dyiiamics  of  z-pinch  plasmas, 
a  description  of  the  transport  of  line  radiation  within  these  optically  thick  plasmas,  and  the 
self-consistent  coupling  of  these  dynamics  to  each  other  and  to  the  fluid  dynamics.  The  model 
employs  the  same  enhanced  transport  coefficients  that  gave  reasonable  agreement  with  the  state  of 
the  plasma  at  the  time  of  stagnation  of  numerous  PI  experimental  argon  and  alumimm  pmches.  ’ 
The  multipliers  on  classical  transport  coefficients  used  in  the  model  are  20  for  resistivity,  30  for 
heat  conduction,  and  40  for  viscosity.  The  calculations  are  begun  with  the  plasma  at  an  average 
radius  of  1.25  cm.  93  %  of  its  mass  is  concentrated  in  a  0.05  cm  sized  shell;  the  remaining  mass  is 
distributed  in  6  inner  cells  inside  the  shell  with  a  mass  density  that  decreases  exponentially  towards 
the  axis.  The  exterior  region  will  be  referred  to  as  the  shell  and  the  interior  region  as  the  core.  This 
initialization  procedure  provides  the  calculations  with  numerical  stability,  and  it  also  simulates  the 
back-pressure  effects  of  plasma  mass  that  is  blown  towards  the  axis  by  shocks,  J  x  B  forces,  or 
gas  jets.  It  is  important  for  modeling  the  core  physics  that  both  ion  and  electron  temperatures  be 
included  in  the  MHD  equations.  If  there  are  no  energy  sources  coupled  into  a  local  plasma  region, 
then  the  ions  and  electrons  would  equflibrate  their  temperatures  typically  on  sub-nanosecond  tune 
scales  in  the  plasma  core;  however,  the  ions  are  preferentially  hea^  rapidly  by  fliermalized  kinetic 
energy  to  hotter  temperature  than  the  electrons.^^  TTiis  nonequilibrium  condition  affects  the  core 

physics  in  that  internal  energy  residing  in  the  ions  cannot  be  radiated  until  it  is  tra^erred  to  the 
electrons.  This  effect  is  pronounced  in  these  studies  because  of  the  way  viscosity  is  enhanced  to 
phenomenologically  model  turbulence. 

In  these  calculations  a  large  parameter  range  in  (t74nass,Z)  space  was  examined.  Thus,  in 
order  to  TniT,inii7/>  the  zone  to  zone  coupling  required  to  transport  the  large  number  of  M-shell, 
L-sheU  and  K-shell  emission  lines  of  the  moderate  Z  elements,  it  was  necessary  to  limit  the  number 
of  zones  used.  On  the  other  hand  one  cannot  do  a  single  zone,  single  temperature,  calculation  by 
averaging  over  all  the  gradients  in  the  plasma  when  in  some  calculations  (and  experiments)  that 
emission  is  produced  only  in  the  hot  core  of  the  plasma,  which  may  contain  less  than  5  %  of  the 
total  mass.  A  credible  calculation  requires  transporting  enough  spectral  lines  to  do  a  reasonable 
job  of  modeling  the  radiative  loss  rates^*  and  at  least  enough  computational  zones  to  reasonably 
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resolve  the  gradients  in  temperature  and  density  that  occur  in  the  plasma  core.  It  is  also  neces^ 
to  resolve  the  outer  boundary  of  the  pinch  because  that  is  where  the  current  flows  and  pene^tw 
into  the  plasma.  Too  course  a  zoning  in  this  region  averages  over  too  much  of  the  colder  mtenor 
plasma  causing  lower  electron  temperatures  to  be  calculated,  which  m  turn,  causes  too  much  fie 
penetration  and  ohmic  heating.  Empirically,  we  found  that  an  adequate  descnpuon  of  the  effect 
of  plasma  gradients  on  K-sheU  emission  could  be  obtained  usmg  15  radial  zones  m  the  alu^^ 
calculation  and  20  zones  in  the  argon,  titanium,  and  krypton  calculations.  Note,  therefore,  for  the 
aluminum  calculations,  the  core  region  contains  1 1  %  of  the  total  plasma  mass. 


per  ion  to  The  total-energy-per-ion  consists  of  the  sum  of  kinetic,  internal  and  pnor  radiated 

energies  of  the  plasma  at  the  time  that  the  current  is  cut-off  divided  by  the  number  of  ions.  As  m  our 
past  work,  these  1-D  calculations  are  used  to  explore  how  effectively  the  thennalization  of 

kinetic  energy  acts  to  drive  the  plasma  on  axis  into  the  K  shell  and  how  effectively  the  to^-enerp 
per-ion  is  converted  into  K-shell  radiation.  With  the  increased  plasma  back-pressure  that  rpults 
from  enhanced  viscosity  and  resistive  heating,  a  large  portion  of  the  energy  coupled  to  the  plpma 
is  in  the  form  of  compressional  work  and  ohmic  heating.  In  order  to  minimize  the  effects  of  ^ 
back-pressure  the  calculations  are  designed  so  that  the  specific  t]*  of  mterest  is  achieved  m  a  0-D 
slug  model  calculation"*  when  the  outer  radius  of  the  plasma  is  a  substantial  distance,  0.5  cm,  from 
the  axis.  In  the  aluminum  calculation,  this  resultsfin  ~  50  %  of  the  plasma  energy  bemg  kmetic 
when  n*  =  1  at  the  time  of  current  termination.  For  all  the  other  higher  Z  or  higher  rj*  calculations, 
the  proportion  of  kinetic  energy  was  higher  than  50  %.  For  the  largest  veloaty  of  krypton 
at  T7*  =  1  the  proportion  was  ~  75  %.  Thus,  calculations  based  on  =  4,  m  earher  work  with 
aluminum  and  argon,*  coixespond  approximately  with  calculations  having  tj*  =  6. 

•  Note,  care  must  be  taken  when  comparing  the  results  of  this  work  with  experiments  or  with 
other  calculations,  not  only  because  of  the  model  limitations  listed  above,  but  also  because  it  is  easy 
to  imagine  pinch  configurations  that  would  not  have  as  large  a  proportion  of  the  total  mput  energy 
being  kinetic  as  is  modeled  here  (such  as  those  having  initial  outer  radu  very  close  to  the  axis  or 
uniformly  distributed  gas  loads).*®  Alternatively,  one  can  imagine  other  pmch  configurpons  that 
achieve  higher  proportions  of  kinetic  energy.  Even  though  it  is  determmedm  this  study  that  the  old 
scaling  laws®  do  not  well  represent  the  low-??  regime,  we  made  no  attempt  to  reformulate  seating 
in  this  regime  because  such  a  model  would  only  be  quantitatively  applicable  to  a  limited  ranp  of 
experiments.  This  point  was  also  made  in  Ref.  (5),  where  &e  K-shell  yield  scalmg  of  alummum 
plasmas  was  determined  for  a  class  of  soft  implosions  (??/??*  =  2/3). 

C,  Results  and  Discussion 

Some  of  the  important  features  of  K-shell  yield  scaling  as  a  function  of  mass  loading  for 
„  =  4  (n*  w  6)  are  mustrated  in  Fig.  1,  which  is  taken  from  the  argon  soft  implosion  calculations 
discussed  in  Ref.  1.  The  mass^  scaling  regime  is  characterized  by  a  <  30  %  conversion  of  coupkd 
energy  into  K-shell  emission.  This  inefficiency  occurs  when  the  radiative  cooting  rates,  w^ch 
scale  as  the  square  of  the  ion  density  (Nf )  or  as  mass^  for  low  mass  plasmas,  are  too  small  to 
effectively  influence  the  dynamics  of  the  pinch.  Under  these  low  mass  conditions,  the  stagnation 
is  characterized  by  an  elastic  collision  whereby  the  plasma  runs  in,  converts  its  kinetic  energy  into 
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internal  energy,  and  then  expands  to  essentially  the  same  kmetic  energy  Ae  ^ 

phase.  The  inability  of  the  plasma  to  radiate  at  maximum  effiaency 

plasma  compression  more  difficult  as  is  characteristic  of  strong  shock  behavior  widi  T  >  ^ 

h  the  ratio  of  specific  heats  at  constant  pressure  to  constant  volume).  As  the  mass  is  mcrease  , 

radiation  cooUng  becomes  sufficiently  large  that  a  signific^t  fraction  of  the 

radiated  away  during  the  stagnation  phase  and  the  time  the  plasma  is  on  axis,  thereby  mar  y 

reducing  the  pinch’s  rebound  dynamics.  K-sheU  emission,  in  this  regime,  scal^  ^c^ly 

input  eLgy,^“  which  is  proportional  to  the  load  mass.  In  this  case,  over  30  %  o^e  kmetic 

energy  coupled  to  the  load  before  stagnation  can  be  emitted  as  K-sheU  radiation. 

which  the  transition  between  mass^  and  mass  scaling  occurs  is  denoted  the  m^break  point 

Because  this  transition  point  is  not  as  clearly  defined  for  “soft”  implosions  as  it  is  for  hard  ,  for 

the  purposes  of  discussion  in  this  study,  the  least  efficient  scaling  regime  is  the  mass  regime  and 

a  more  efficient  regime  begins  when  it  is  obvious  that  the  paling  has 

TTierefore,  the  criterion  of  achieving  >  30%  conversion  of  kmetic  energy  mto  K-sheU  photom  ^ 
being  a  condition  of  reaching  the  efficient  scaling  regime,  which  was  used  m  earlier  work, 
not  strictly  adhered  to  in  this  study. 

It  is  easier  to  insure  efficient  production  of  K-shell  emission  by  imploding  a  plasma  load  in 
the  high-mass  high-t?  scaling  regime;  however,  except  for  low-Z  loads,  the  “7^ 

required  to  operate  in  this  regime  quickly  becomes  prohibitive  as  Z  mcreases.  This  is  illustrated  m 
Fk  2  where  calculated  load  energies  required  to  obtain  10  %  conversion  of  coupled  ener^  mto 
K-shell  X  rays  are  shown  as  a  function  of  Z.  For  titanium  it  is  estimated  that  approi^ately  100 
kJ/cm  of  energy  must  be  delivered  to  the  load  while  krypton  requires  7  MJ/cm.  For  dements  iwth 
higher  atomic  number  than  Z  =  22  (titanium),  more  energy  is  required  to  obtam  10%  efficiency  than 
cm  be  deUvered  by  existing  pulse  power  machines.  Because  of  this  ener^  Citation  modeiate-Z 
experiments  performed  on  existing  pulse  power  generators  will  take  place  m  low-iy  scalmg  regimes. 

K-shell  aluminum  yields  as  a  function  of  load  mass  for  various  t?*  are  displayed  in  Fig.  3. 
Each  soUd  fine  curve  corresponds  to  a  fixed  value  of  tj*  attained  at  the  time  of  current  termination. 
The  solid  curves  in  this  figure  shows  that  an  ineffident  scaling  regime  extends  to  much  larger  mass 
loads  than  predicted  by  earlier  modeling  worL  For  example,  the  mass  scalmg  regime  exten^ 
to  between  100  and  200  /ig/cm  for  all  values  of  v*  between  1  and  fi^^lier  predretions  fm  e 
Mfcp’s  are  10, 35,  and  124  fig/cm,  for  ,  =  1  (2),  2  (3),  and  4  (6).  •  ^  n^bers  m  Par®“7?7a 
are  roughly  the  equivalent  v*  values.  The  earlier  predictions  are  based  on  hard  implosion  sc^g 
laws  that  were  modified  by  a  factor  of  six  more  mass  to  match  soft  implosion  17  -  4  calcul^ons^ 
Obviously,  our  present  calculations  predict  that  more  than  this  f^tor  of  six 
implosion  scaling  relations  is  necessary  in  order  to  capture  the  dynamics  of  K-shell  yield  scalmg 

in  the  low-»7*  regime. 

The  dashed  curves  in  Fig.  3  show  the  K-shell  yield  as  a  function  of  mass,  but  ener^ 

deUveredtotheloadasaparameterratherthanr,-.  These  curves  have  a  consistent  trend.  T^^show 

thatthepeakinyieldproductionriseswithiy*  as  theenergy  delivered  to  Aeloadmcrea^.  Tlu^,  at 

small  energy  coupling,  the  optimal  choice  is  to  take  advantage  of  the  JST-  dependence  of  K-shell  Im 
radiation  id  to  Lplode  as  much  mass  as  possible  with  17*  «  1-  However,  as  mass  mcre^s, 
two  mechanisms  begin  to  limit  K-sheU  emission  from  r,*  =  1  imploaons.  First,  temperatoes 
fall  because  the  plasma’s  ability  to  radiate  rises  more  rapidly  than  its  ability  to  thermalize  kmetic 


rt  K  chpll  line  ooacitv  (including  that  due  to  innershell  absorption)  also  rises. 

of  fj*). 

The  dvnamics  of  (he  first  meohanism  are  fflostrated  in  Fig.  4.  which  compares  one  power 
inonm  w^er  outputs  as  a  function  of  load  mass  for  ,•  =  1:  the  peak  rate  of  ^ 

^  total  radiative  power,  and  the  peak 

mass  scaling  regime  both  radiative  powers  increase  with  mass  ;  where^,  the  rate  at  w 

oftheplasma««co™^^^ 

Cerra“-uli:se;ntodropasa~fm^^^^^^^^^^ 

mass  represents  the  source  region  for  90  %  of  the  K-sheU  rauianom  oronortion  of 

time  integrated  emission  than  absorption  are  considered  as  part  of  this  gi  ^  ^  -i  ui,.  tn 

te  ntema^mllg  in  the  K-sheU  increases  with  ,•  because  of  the  addmond  energy  ^vad^  e  to 

wSSSmtheKsheU.  Fig.  ^  ^Ilsb 

fractional  mass  emitting  in  the  K-shell  decreases  with  mass.  The  source  of  the  K-sheU  emi^i 
Se  ltg“  mass  plasmas  is  the  core  and  nearby  regions.  The  K-shell  spectmna  transtnons  from 
Lg  P^lminanSy  lines  for  the  low  mass  implosions  to  being  mosfly  recombmatton  radtauon 

for  the  high  mass  implosions. 

For  large  mass  loads,  the  amount  of  K-shell  radiation  that  ultimately  es^pes  P" 
much  less  than  that  escaping  the  source  region.  This  reducnon  occurs  because  the  K-sheU  phmom 
generated  in  the  hot  interior  source  region  are  absorted  by  innemheU  ^ 

:rrp“»^"y  50  »  when  innersheU  photoahsorption  was  taken  to  a^o»n 
r  1  ^  1  sac  tViam  m*  =  1  the  cffcct  of  innersheU  absorptiou  IS  not  as  pronounced  for  any 

spedSmals  load  becaL  there  is  less  cold  absorbing  plasma  (Fig.  6).  ^°’‘^SherZ  element  to 

specmc  mass  loau  emission  but  now  at  mass  loads  that  scale 

^™otoraLiinum  mass  where  diis  effect  occurs.  This  finding  reflects  the  Z  dependence 
of  the  innershell  absorption  coefficient,  i.e., 

&  =  {glcm?), 

wherckisaconstanhZistheatomicnumherofd.eabsorhertoAisd.ephot««w^^^^^^ 

i  ^  7-2  for  above  K-edge  threshold  photons,  the  Z  dependence  of  o  is  Z  .  These  mommtm 
Lcufations  suggest  that  the  effects  of  innersheU  absorption  on  K-sheU  mdtanve  yield  can 
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minimized  by  using  a  mixture  plasma  that  consists  of  a  different  Z  extenor  plasma  whose  K-edge 
Ues  far  from  the  K-shell  photon  energies  of  the  interior  plasma.  Recently  performed  experiments,^^ 
which  were  based  on  the  idea  of  using  plasma  mixtures  to  reduce  the  overall  opacity  of±c  pl^ma, 
provide  evidence  that  the  effects  of  reabsorption  can  be  large  and  can  be  lessened  by  this  techmque. 
As  Fig.  5  shows,  radiation  generated  by  the  interior  plasma  becomes  trapped  by  the  extenor  plasma 
and  causes  the  shell  temperature  to  rise  slowly  with  mass  when  the  opacity  is  suffiaently  large  m 
large  mass  load. 

The  yield  predictions  of  Fig.  3  are  normalized  to  the  coupled  energy  to  show  the  efficiency 
of  converting  coupled  energy  into  K-shell  emission  as  a  function  of  mass  load  m  Fig.  »•  Sjnee 
ti*  =  1  implosions  barely  have  enough  energy  to  ionize  the  entire  plasma  mto  the  K-shell  (and 
this  consideration  neglects  radiation  losses)  one  might  expect  a  more  precipitous  drop  m  the  peak 
efficiency  of  producing  K-shell  emission  as  r/*  approaches  1.  The  reason  this  drop  does  not  occur 
is  that  the  geometry  and  implosion  dynamics  of  the  z-pinch  focuses  the  implosion  energy  towards 
the  axis  so  that  the  i/*  values  in  the  core  of  the  pinch,  which  is  where  the  K-shell  photons  are 
generated,  can  be  much  greater  than  one. 

It  was  mentioned  earlier  that  the  inability  of  a  plasma  to  radiate  fast  enough  relative  to  the  rate 
kinetic  energy  is  thermalized  leads  to  an  implosion  dynamics  that  is  characteristic  of  strong  shock 
behavior,  consequently,  it  is  difficult  to  compress  a  low  mass  plasma.  This  mcompressibility  is 
illustrated  by  Fig.  9.  In  this  figure,  the  mass  and  time  averaged  shell  and  core  ion  densities  and 
electron  temperatures  are  shown  as  a  function  of  17*  for  a  100  fig/cm  aluminum  load,  '^e  ^e 
average  is  carried  out  over  the  full  width  half  maxmum  of  the  K-shell  emission  pulse  width.  Fig. 

9  shows  that  the  core  ion  density  is  insensitive  to  the  velocity  of  the  implosion,  and  to  a  les^r 
extent  the  same  is  true  for  the  shell  ion  density.  If  the  implosions  are  modeled  using  a  circuit  model 
appropriate  for  a  real  machine,  and  the  current  is  not  terminated  prior  to  stagnation,  one  finds  that 
the  energy  that  can  be  coupled  into  a  low  mass  load  is  significantly  limited  by  the  mability  of  the 
current  to  do  compressional  work  on  the  plasma. 

Figs  10  and  1 1  display  K-shell  yields  as  a  function  of  mass  load  and  rj*  for  argon  and  titamum. 
In  these,  as  in  the  aluminum  calculations  (Fig.  3),  1)  there  is  a  characteristic  mass^  scaling  regime 
that  traiisitions  to  a  more  efficient  scaling  regime  as  the  mass  is  increased  and  the  radiation  rates 
begin  to  compare  with  the  rate  kinetic  energy  is  thermalized,  2)  the  transition  from  an  mefficient 
to  an  efficient  scaling  regime  occurs  at  much  larger  mass  loads  than  predicted  by  earlier  scalmg 
laws  3)  yield  predictions  for  a  fixed  amount  of  coupled  machine  energy  (dashed  curves)  show  that 
the  optimal  load  choice  for  producing  K-sheU  emission  favors  larger  masses  and  lower  implosion 
velocities  Gower  17*)  at  low  machine  energy  coupling,  and  4)  the  peak  efficiency  of  produ^g 
K-shell  emission  diminishes  with  rj*  (not  shown).  The  peak  efficiencies  for  argon  with »?  —  ’  • 
3,  and  4  are  respectively,  24, 27, 32,  and  33  %.  For  titamum  with  17*  =  1, 2,  and  3  the  efficiencies 

are  18, 23,  and  27  %. 

The  cross-hatched  regions  shown  in  Fig.  11  represent  predicted  operating  regions  for 
optimizing  K-shell  X-ray  production  on  Double  EAGLE,  Blackjack  5,  Phoenix,  and  Saturn 
class  machines.  They  represent  good  regions  in  which  the  Z-scaling  of  K-shell  emission  can 
be  investigated  experimentally.  Titanium  is  a  good  load  choice  because:  1)  its  atomic  number 
Z=22  is  significantly  above  that  of  argon,  which  has  been  extensively  explored,  2)  a  large  mass  ot 
titani^  can  be  ionized  into  the  K-shell  by  existing  pulse  power  machines,  and  3)  titanium  can  be 
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“n  a^s  lL  1-D  symmetry.  Recent  Al,  and  A1  coated  with  Mg  experiments  perform^ 

return  geometry  by  putting  the  return  current  path  at  a  larger  radius. 

H/lto  and  H-like  ioniaation  states  in  a  time  short  enough  to  come  mto  CRE  at  the  calcul^ 

He-ii^  ana  n  —  1  ralculations  these  conducive  conditions  exist  for  a 

(TArp  rpfnon  of  the  low  mass  alummuni  loads  (ri§*  to  y  *^4.^  tr 

L  notTng  enoSi  for  the  lower  mass  krypton  loads  (Fig.  12).  The  penod  for  to^g  mm  the  K 

shell  is  largely  influenced  by  the  total  radiation  rate  and  the  rate  kinetic  energy  IS  erm 

For  a  fixed  electron  density,  this  time  increases  with  the  atomic  number  of  the  load  maten^_ 
elect^nimp^^^ 

time  reauired  to  reach  the  next  sequential  iomzation  stage  is  given  by  r  -  {N.  X )  ,  w  e 

is  dm  decnon  density  and  X  is  the  collimonal  dre 

TLi-He  is  written 

F.,(^)N.y'  see,  (2) 


r  >  TLi-He 


=  ( 


XI 


,1/2 


rSetr^  coupled  m  the  load)  can  he  esdmated  by  «  W 


10^®cm  ^ 


an  lomzanon  uiuc  ui  j.  »e-o  - - - r-  -  .  ,  _. 

for  ainminun.  and  approximately  25  times  this  density  for  krypton. 
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Besides  the  necessity  of  higher  electron  density,  more  sequential  lom^hon,  and  short 
ionization  times,  the  enhanced  viscosity  employed  in  our  model  makes  it  more  difficult  m  theo^ 
for  the  higher  Z  elements  to  ionize  into  the  K  shell.  This  is  because  the  enhanced  ^scosity  produces 
core  ion  temperatures  during  the  stagnation  phase  of  the  implosion  that  are  orders  ® 

larger  than  the  electron  temperatures.  Significant  energy  residmg  m  the  ions  makes  it  “ore  difficult 
to  compress  the  plasma  which  behaves  like  an  ideal  gas  until  the  enerp  is  tr^ferred  to  the 
electro^  where  it  can  be  radiated  out  of  the  plasma.  Lower  compression  imphes  JV*  whic 
in  turn  implies  longer  ionization  times  and  also  longer  time  for  the  electrons  to  equili  ra  wi 
ions  since  the  ion  electron  heating  rate  scales  as  Ng . 

D.  Summary  and  Conclusions 

Recent  investigations  into  phenomenological  turbulence  modeling  showed  ^ 

enhancement  of  viscosity,  heat  conduction,  and  resistivity  in  1-D  MHD  calculations  couW  lead 
to  reasonable  agreement  with  the  average  stagnation  state  of  the  experimental  plasmas.  ’ 
particular  stagnation  ion  densities  were  reduced  by  several  orders  of  magmtode  over  the  de^i^es 
predicted  using  classical  transport  coefficients.  This  same  modeling  was  used  m  the  workpresented 
here  to  further  investigate  the  Z-scaling  of  K-shell  emission  from  moderate-  Z  plasmas  at  low-i?. 
The  low-77  regime  is  investigated  because  of  limits  on  energy  available  from  current  pulse  pow^ 
machines  to  drive  loads  with  higher  atomic  number  than  Z=22  (tit^um) 

study  have  identified  several  important  effects  that  will  influence  the  sc^g  of  K-sheU  emission 
on  these  machines.  First,  low-17  plasmas  remain  in  an  inefficient  mass  scaling  regime  at  much 
higher  masses  than  predicted  by  our  earUer  work.  Secondly,  opacity  and  innersheU  absorption 
limit  the  K-sheU  radiative  capabilities  of  the  higher  mass  loads.  Thirdly,  there  is  a  Z  dependent 
lag  in  ionization  into  the  K-shell  that  increases  with  the  Z  of  the  plasma  and  adversely  affects  the 

K-shell  emission. 

It  was  shown  in  Ref .  1  that  there  is  a  factor  of  6  upward  shift  in  the  minhnunt  load  needed 

for  producing  efficient  K-sheU  yields  in .,  =  4  “soft”  as  opposed  to  W  implosion  calcu^o^ 
HoLer,  one  would  mistakenly  ptedictMt,(,)  “soft”  =  6x  Ms,(l)  for  low-,  implosto^ 
This  erroneous  result  would  also  assume  that  one  can  extrapolate  the  hard  unplosion  sc^g 
Mr  (n1  to  low  It.  We  found  that  1-D  soft  implosion  calculations  do  not  support  tins  assumpnon. 
We  fomid  that  low-,  plasmas  remain  in  an  inefficient  scaling  regime  at  much  higher  masMS  than 
r^Ld  earUer.  e.g  M.p(,  =  1)  “soft”  »  60x  Ms,(,  =  1)  “hard”  for  alummum.  Beca^ 
L  calculations  predict  that  mote  mass  will  he  needed  to  obtain  an  efficient  generanon  of  K-shell 
emission  than  earher  predictions  they  also  predict  more  energy  will  be  required. 

Our  calculations  consistently  show  that  the  optimal  conditions  for  generating  K-shell  photons, 
for  a  specific  load  energy  coupling,  transition  from  low- 17*  higher  mass  unplosion  conditions  to 
higher^*  lower  mass  conditions  as  the  coupled  energy  is  increased.  TWo  mechamsms  that 
tli  transition  are:  1)  bulk  plasma  temperatures  falling  with  mcreased  mass  as  the  radiative  coo  g 
rate  (scales  with  maaa^)  begins  to  compare  with  the  rate  totic  energy  is 
with  maas),  and  2)  K-shell  photons  initiaUy  emitted  from  the  core  are  mcreasmgly  absorbed  by 
innersheU  processes  in  colder  outer  regions  of  the  plasma  as  the  mass  is  mcreased. 

In  both  “soft”  and  “hard”  implosion  calculations  the  core  region  is  much  hotter  than  the 
surrounding  sheU  plasma  during  the  stagnation  phase  of  the  implosion.  This  hot  core  is  more 
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extensive  and  its  ion  temperatures  am  substantially  higher  than  the 

••soft”  caloulatioos  are  performed.  Hie  Mg  coated  A1  wne  expemnents 

_ . ,  _ : _ n„i  tVitit  the  statp,  of  the  core  olasma  is  much  different  than 


rr^dtS"  i.  was  determined  that  the  core  region  was  approximately 

twice  as  hot  as  the  shell.  In  order  to  gain  better  knowledge  of  how  K-sheU  radiabon  is 

is  clear  that  more  experiments  are  needed  to  observe  and  to  quantify  the  state  of  the  core  and  sheU 

regions. 

If,  for  example,  on  a  specific  pulse  power  machine  it  were  possible  to  couple  ^  “"^h  ^ 
energy  into  shorter  length  more  massive  loads  as  into  longer  length  less  massive  loads  and  if 
mX  Ire  energy  Lited  and  forced  to  operate  in  an  inefficient 
element,  then  it  would  be  advantageous  to  implode  a  more  massive  shorter  length  load 
m  ^  advlmge  of  scaling  For  higher  Z  loads  like  krypton,  there  is  a  double  advantage 

for  using  the  largest  mass  per  unit  length.  For  this  element,  yield  scalmg  m 

is  evL  steeper  Im  mass^  (Fig.  12)  because  of  the  of  the  longer  tunes  required  for  lomzation 

into  the  K  sb^U.  For  simUar  reasons,  load  designs  that  promote  higher  mass  and  energy  density  at 

stagnation  such  as  convergent  geometry  loads^^-^^  or  loads  that  implode  onto  a  wire  on  axis,  may 

alsfimprove  the  efficiency  of  generating  K-shell  emission.  In  addition,  improvi^  to 

the  implosion  by  going  to  more  wires  and  more  symmetric  return  current 

plasma  shearing  on  axis  and  the  need  to  enhance  transport  coefficients  m  order  to  at^  I  D  m 

agreement  with  present  day  experimental  stagnation  conditions.  If  this  is  &e  c^.  then  the  ener^ 

riuirements  predicted  here  for  attaining  efficient  K-sheU  emission  may  be  substa^aUy  reduced 

S  become  more  in  accord  with  the  original  scaling  laws.*  These  matters  need  to  be  mvestigated 

experimentally. 

References 

1.  J.  W.  ThomhiU,  K.  G.  Whimey,  C.  Deeney  and  P.  D.  LePell,  Phys.  Plasmas  1, 321  (1994). 

2.  C.  Deeney,  P.  D.  LePell,  B.  H.  Failor,  J.  .  S.  Meachum,  S.  Wong,  J.  W.  Thornhill,  K.  G. 

Whitney,  and  M.  C.  Coulter,  J.  Appl.  Phys.,  75, 2781  (1994). 

3.  C.  Deeney,  T.  Nash,  R.  R.  Prasad,  L.  Warren,  J.  W.  Thornhill,  and  M.  C.  Coulter,  Phys.  Rev. 

A,  44, 6742  (1991). 

4  K.  G.  Whitney,  J.  W.  ThomhiU,  J.  L.  GiuUani,  Jr.,  J.  Davis,  L.  K  ^es,  E.  E. 

Kenyon,  W.  L  Speicer,  J.  A.  Draper,  C.  R.  Parsons,  P.  Dang,  R.  B.  Spielm^,  Nash.  L 

S.  McGum,  L.  E.  Ruggles,  C.  Deeney,  R.  R.  Prasad,  and  L.  Warren,  Phys.  Rev.  E,  50, 2166 

(1994). 

5.  K.  G.  Whimey,  J.  W.  ThomhiU,  J.  P.  Apmzese,  J.  Davis,  C.  Deeney,  P.  D.  LePeU  and  B.  H. 

Failor,  Phys.  Plasmas  2, 2590  (1995). 

6.  K.  G.  Whimey,  J.  W.  ThomhiU,  J.  P.  Apmzese  and  J.  Davis,  J.  Appl.  Phys.  67, 1725  (1990). 

7.  J.  W.  ThomhiU,  K.  G.  Whimey  and  J.  Davis.  J.  (^uant  Spectrosc.  Radiat  Transfer  44,  251 

(1990). 

8.  D.  Duston,  R.  W.  Clark,  J.  Davis,  and  J.  P.  Apmzese,  Phys.  Rev.  A  27, 1441  (1983). 


28 


9.  J.  P.  Apruzcsc,  J.  Davis,  D.  Duston,  and  R-  W.  aark,  Phys.  Rev.  A  29, 246  (1984). 

10.  J.  P.  Apruzese,  J.  Davis,  D.  Duston,  and  K.  G.  Whitney,  J.  C^iant  Spectrosc.  Radiat  Transfer 
23, 479  (1980). 

11.  C.  Deeney,  P.  D.  LePell,  B.  H.  Failor,  S.  L.  Wong,  J.  P.  y^ruzese,  K.  G.  Whitney,  J.  W. 
Thornhill,  J.  Davis,  E.  Yadlowslq^,  R.  C.  Hazelton,  J.  J.  Moschella,  T.  Nash  and  N.  Loter,  Phys. 
Rev.  £,51,4823(1995). 

12.  J.  P.  Apruzese,  R.  W.  Clark,  and  J.  W.  Tliomhill,  J.  Comp.  Phys.  119, 156  (1994). 

13.  J.  Davis,  J.  L.  Giuliani,  Jr.,  and  M.  Mulbrandon,  Phys.  Plasmas  2, 1766  (1995). 

14.  J.  P.  Apruzese  and  J.  Davis,  Naval  Research  Laboratory  Memorandum  Report  No.  5406  (1984). 

15.  R  L.  Cochran,  J.  Davis,  and  A.  L.  Velikovich,  Phys.  Plasmas  2, 2765  (1995). 

16.  This  concept  was  suggested  by  R.  B.  Spielman  of  Sandia  National  Laboratories  in  a  private 
communication. 

17.  B.  CuUity,  in  Elements  of  X-Ray  Diffraction,  (Addison-Wesley,  1978). 

18.  J.  P.  Apruzese  and  J.  Davis,  J.  Appl.  Physics  57, 4349  (1985). 

19.  TTie  fact  that  modeling  time  dependent  ionization  dynamics  is  important  for  krypton  z- 
pinch  implosions  was  suggested  by  T.  NasL  of  Sandia  National  Laboratories  in  a  private 
communication. 

20.  R.  B.  Spielman,  D.  L.  Hanson,  M.  A.  Palmer,  M.  K.  Maatzen,  T.  W.  Hussey,  and  J.  M.  Peek, 
J.  Appl.  Phys.  57, 830  (1985). 

21.  W.  Lotz,  Z.  Phys.  216, 241  (1968). 

22.  J.  H.  Degnan,  R  M.  Lehr,  J.  D.  Beason,  G.  P.  Baca,  D.  C.  Bell,  A.  L.  Chesley,  D.  Dietz,  D.  B. 
Dunlap,  S.  E.  Englert,  T.  J.  Englert,  D.  G.  Gale,  J.  D.  Graham,  C.  D.  Holmberg,  T.  W.  Hussey, 
R.  A.  Lewis,  C.  A.  Outten,  R.  E.  Peterkin,  Jr.,  D.  W.  Price,  N.  R  Roderick,  E.  L.  Ruden,  U. 
Shumlak,  G.  A.  Smith,  and  P.  J.  Tlirchi,  Phys.  Rev.  Lett.  74, 98  (1994). 

23.  I.  V.  Sokolov,  Sov.  Phys.  Usp.  33, 960  (1990). 


29 


kinetic  ener! 


•meats  for  obtaining  10  pete 
a  function  of  atomic  number. 


33 


Figure  So  Ma^  and  time  averaged  clectroE 

regions  of  an  alraniffliffli  tj*  =  1  plasma  as^  a  to 
fuH-wdth-half-masimum  of  the  K-shell  cmisami 


I  tOTpcratuies  and  ion  densities  of  the  core  and  s 
a  fimction  of  mass  load.  The  averaging  is  over 


•I 


T-rtor?  T<^ 


Mass  (|J.g/cm) 


Figure  6.  Fractional  mass  anitting  K-shell  radiation  as  a  function  of  mass  load  and  ri*  for 
aluminum  z-pindi  implosions. 
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in.  nonthermal  radiation  from  imploding  z-pinch  plasmas 

During  the  implosion  and  stagnation  of  a  Z-pinch  plasma  a  fr^tion  of  the  current  is  carried 
by  nonthermal  electrons  that  produce  energetic  x-rays  during  collisions.  Here  we  estinmte  e 
potential  for  producing  nonthermal  bremsstrahlung  below  100  keV  ba^d  on  some  extreme  y 
optimistic  assumptions:  1)  100  %  of  the  current  is  carried  by  nontermal  electrons,  2)  all  of  the 
energy  coupled  to  the  pinch  goes  to  generating  this  current,  3)  The  f^t  electron  energy  is  comple  y 
absorbed  in  the  plasma  (thick  target  bremsstrahlung),  and  4)  there  is  no  opacity. 


The  energy  per  incident  electron  that  is  emitted  as  bremsstrahlung,  Q,  is  given  by 


dE 


di 


■57  Rad  -r  Jo 


where  the  kinetic  energy  of  the  incident  electron  is  E  (See  Evans  -  “The  Atomic  Nucleus”). 


^\B^d=N 

0 

.t 

is  the  relativistic  radiative  loss  rate  per  unit  path  length  ds,N  is  ion  density,  and  xH  is  the 
radiative  cross  section  (see  Jackson  -  “Classical  Electrodynamics”,  1962  ed.). 


is  the  relativistic  ionization  loss  rate  per  unit  path  length  ds  (see  Dalgamo  -  “Atomic  and  Molecular 
Processes”). 


The  radiative  yield  in  a  frequency  range  which  is  due  to  a  single  incident  electron  with 

energy  E,  is  calculated  from 

J 


F(«ui,<«;2)  =  - 


The  percentage  of  the  incident  electron  energy  that  appears  as 
range  wi  —  W2  is  then 


y’(wi,w2)  _  _ 

E  ^  E  ^ 


bremsstrahlung  in  the  frequency 
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Assuming  that  all  of  the  machine  energy  that  is  coupled  into  the  pinch  goes  into  generating 
nonthermal  electrons  at  a  given  incident  energy  E,  then  the  total  radiative  yield  in  the  frequency 
range  (wi  —  wj)  normalized  to  the  machine  energy  is  also  given  by 


Total  Radiation{u}i ,  to>2  ) 
Machine  Energy 


This  is  the  best  that  can  be  done  because  it  also  assumes;  1)  100  percent  of  the  current  is  earned  by 
nonthermals,  2)  any  arbitrarily  large  voltage  V  or  reflexing  that  is  needed  to  produce  the  electrons 
with  energy  E  can  occur  across  the  pinch,  i.e.,  it  allows  for  voltage  or  power  amplification,  and  3) 
no  radiative  opacity.  The  results  of  this  work  as  appUed  to  a  general  pulse  power  machine  and  a 
tungsten  target  is  shown  in  Fig.  1  and  results  specific  to  Double  EAGLE  and  Saturn  are  shown  m 
Figs.  2  and  3,  respectively. 

Nonthermal  Electron  Generated  K-alpha  Emission 


We  applied  the  same  analysis  as  used  for  thick  target  bremsstrahlung  to  K-alpha  emission  from  a 


tungsten  plasma. 


where. 


\x-a'=  Na’{E,Z)  X  ionization  potential 


tr  is  the  cross  section  for  knocking  out  an  innershell  electron  (Coulomb-Bom-exchange  calculation), 
which  is  dependent  upon  the  atomic  number  Z  of  the  absorber  and  the  incident  energy  E  of  tiie 
electron.  N  is  the  ion  number  density.  The  K  fluorescence  yield  is  taken  as  unity  for  tungsten  and 
the  ionization  potential  is  69.5  keV.  The  results  for  nonthermal  electron  generated  K -a  emission 
from  a  general  pulse  power  machine  are  shown  in  Fig.  4. 


We  conclude  that  our  (overly  optimistic)  estimate  shows  that  there  is  a  potential  for  2  % 
conversion  of  machine  energy  into  nonthermal  electron  generated  bremsstrahlung  below  100  keV 
from  tungsten  (conversion  efficiency  is  proportional  to  Z).  In  reality  this  estimate  is  probably 
generous  by  over  an  order  of  magnitude.  This  conclusion  is  supported  by  (Tiger  code)  three 
dimensional  Monte-Carlo  electron,  photon  transport  calculations  (thick  target)  that  show  conversion 
efficiencies  that  are  two  orders  of  magnitude  below  our  simple  optimistic  estimates.  Based  on  this 
analysis  it  does  not  appear  that  nonthermals  are  an  interesting  source  for  producing  radiation  in  the 
20  - 100  keV  range. 
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K-alpha  Radiation  /  Machine  Energy 


K-alpha  Results 

Under  the  assumptions  of  optically  thin  opacity  and  100  % 
conversion  of  machine  energy  to  nonthermal  electrons. 


wire  by  the  wsnnuer  argon  plasiiisu 


me _ 

below  iiiTable  1. 


Puff  on  Wire 


Energy  coupled  to  load  (kJ/cni) 

15 

20 

Total  radiative  yield  (kJ/cm) 

7 

13 

Peak  radiative  power  (TW/cm) 

0.3 

0.5 

Radiation  >  3  keV  (kJ/cm) 

0.5 

0.9 

Peak  >  3  keV  power  (TW/cm) 

0.04 

0.08 

Peak  core  <Ni>  (cm"*) 

2x10^® 

4x10^® 
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Peak  core  <Te>  (keV) 
Peak  core  <Ti>  (keV) 


1.8 


1.6 

20 


40 


Table  1. 


where  <>  represents  mass  and  time  averaged  quantities  over  the  plasma  that  is  emitting  m  e 
K-shell  ~  10  %  of  the  plasma  (located  in  the  core).  These  results  clearly  show  that  the  PO 
configuration  for  the  chosen  conditions  couples  more  energy  into  the  load,  radiates  almost  twice  ^ 
much  in  the  total  yield  and  >  3  keV  regime,  and  is  twice  as  dense  in  the  core  than  a  pure  gas  putt. 


An  experimental/theoretical  investigation  into  the  puff  on  wire  configuration  was  also  miti^d 
with  PI  The  experiments  were  performed  using  an  argon  gas  puff  and  titamum  wes^ 
experimental  pinches  were  4  cm  in  length  and  the  nozzles  were  1.25  cm  in  radius.  The  Double 
EAGLE  circuit  model  consisted  of  0.3  ohms  in  series  with  20  nh  inductee  m  senes  If 

The  open  circuit  voltage  profile  is  described  by  Voc  =  2.35MV’ (1  +  (t/320)  )am(7rt/160),  where 


i  is  inns. 


Table  2.  shows  a  summary  of  results  from  these  calculations.  There  are  results  for  puff  only 
calculations  of  100, 150, 200  and  250  ng/cm  as  weU  as  results  for  these  same  gas  puffs  impmging 
upon  0.8  mil  and  1.0  mil  titanium  wires.  Note,  we  originally  thought  the  pinches  were  3  cm  m 
length  so  most  of  the  calculations  were  done  at  this  length.  We  have  not  yet  gone  back  and  repeated 
the  calculations  at  4  cm.  However,  we  did  perform  h  calculation  for  a  4  cm  long,  150  ^p/cm  argon 
puff  onto  a  0.8  mil  titanium  wire  load  for  which  it  appears  to  first  order  fiiat  die  radmtive  energy 
that  is  lost  per-unit-length  is  made  up  by  having  a  longer  pinch.  Also  the  implosion  time  was  ns 
longer  for  the  4  cm  length  load  because  of  its  larger  inductance. 

Fig  2  shows  a  comparison  between  experimental  and  calculated  >  3  keV  yields  for  puff-only 
and  puff  onto  a  0.8  mil  titanium  wire.  The  implosion  times  are  all  fiduciated  with  the  linear  rise  of 
the  current  and  4  ns  has  been  added  to  the  calculated  implosion  times  in  order  to  correspond  more 
closely  with  the  experimental  times.  From  this  figure  we  see  that  there  is  reasonable  agreement 

betweencalculationsandexperimentsforthepuff-onlyconfiguration.However,moreexpenme^ 

data  is  needed  at  both  longer  (more  massive)  and  shorter  (less  massive)  implosions  times.  The 
experimental  and  theoretical  yields  are  also  in  reasonable  accord  when  the  0.8  mil  wire  is  on  axis 
The  ACE  4  calculation  that  was  discussed  above  and  the  lower  implosion  time  results  (Fig.  ) 
show  that  at  small  mass  loads  the  POW  configuration  out  performs  the  puff-only  configuration. 
Unfortunately,  there  were  no  puff-only  experiments  performed  at  low  imptosion  times  to  ver^ 
these  calculated  trends.  This  is  an  important  trend  that  should  be  mvestigated  expenmentaUy 
because  the  POW  may  be  a  technique  for  obtaining  maximal  K-shell  ^fields  from  elements  that  c^ 
only  be  imploded  in  the  inefficient  scaling  regime  because  of  machine  energy  linutations,  i.e.,it 
may  improve  maximal  titanium  yields  on  Double  EAGLE  and  Saturn,  krypton  yields  on  a  Jupiter 
class  machine,  or  argon  yields  on  ACE  4. 


Fig.  3  shows  a  comparison  between  experimental  and  calculated  >  3  keV  yields  for  puff-only 
and  puff  onto  a  1.0  mil  titanium  wire.  Because  of  the  large  differences  between  calculated  and 
experimental  yields  at  the  larger  implosion  times,  one  has  to  question  the  validity  of  way  the 
interaction  of  the  inner  puff  gas  and  titanium  wire  is  bemg  modeled.  One  reason  the  calculated 
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yields  drop  off  at  the  larger  implosion  times  is  that  the  1  mil  titanium  ^ 
i:g„nplasLastotheeamnt^.ti.e.^ 

^e“  !n  ti^  is  possibly  an  indication  that  we  are  over«tima^g  tim 

into  the  titanium  wire  or  that  the  experimental  core  densities  are  not  as  l^e  as  p  edtcted.  A  closer 

examination  of  the  experimental  data  is  needed  m  order  to  better  asses  the  mo  e  . 

In  addition  to  optimizing  K-shell  emission  for  low  mass  loads,  improt^g  simulation  fideUty 
was  another  reason  for  investigating  the  POW  configuration.  ^  goal  o  e 
improve  the  yield  above  4  keV.  Ihe  data  are  still  being  analyzed  but  one  encoura^g 
Zlt  is  thatl  best  puff-only  experiments  produced  100  J  of 

whereas  the  91  ns  implosion  time  puff  onto  a  0.8  mU  titanium  wire  load  produced  m  excess  of  1 
kJ  of  radiation  in  this  range.  The  yield  above  3  keV  may 

nuff-onlv  configuration  for  this  load  -  needs  further  experimental  mvestigation.  O^e  c^culations 
performed  the  one  that  best  represents  this  experiment^  load  is  a  3  cm  length  1 

Duff  onto  a  0.8  mil  titanium  wire  calculation.  The  calculated  yield  m  the  4.5  -  5.0  keV 

L  from  this  load  and  it  is  225  J  from  the  puff-only  configuration.  This  result  also  indicates  Aat  4 
core  ion  densities  may  be  overpredicted  since  the  calculated  puff-only  yields  m  the  range  4.5-5.0 
keV  are  larger  than  measured. 


50 


PI  Argon  Puff  onto  Ti  Wire  Shot  Matrix 

Machine  —  Double  EAGLE 
Charging  voltage  60  kV 

1.25  cm  radius  nozzle  Ar„nn  only 

3  cm  length  pinch 


>3keVYield  (kJ/cm) 

>  4  keV  Yield  (kJ/cm) 
Implosion  time  (ns) 

Coupled  energy  (kJ/cm) 
* 
n 

Peak  kinetic  e.  (kJ/cm) 

n 


100  mg/cm 
3.0 
0.4 
102 
24 
2.5 
8.4 
0.9 


150  mg/cm 

5.0 

0.7 

no 

29 

2.0 

10 

0.7 


200  mg/cm 

6.0 

0.9 

118 

34 

1.9 

11 

0.6 


1  mil  Ti  wire 


>  3  keV  Yield  (kJ/cm) 

>  4  keV  Yield  (kJ/cm) 
Implosion  time  (ns) 

Coupled  energy  (kJ/cm) 
* 

Peak  kinetic  e.  (kJ/cm) 

•n 


0.8  mil  Ti  wire 


>  3  keV  Yield  (kJ/cm) 

>  4  keV  Yield  (kJ/cm) 
Implosion  time  (ns) 

Coupled  energy  (kJ/cm) 

* 

T] 

Peak  kinetic  e.  (kJ/cm) 

n 


250  mg/cm 

6.0 

0.9 

125 

36 

1.6 

11.5 

0.5 


0.4 

0.02 

125 

36 

1.6 

11.5 

0.5 


3.0 

1.0 

125 

36 

1.6 

11.5 

0.5 


Note  quoted  implosion  times  are  fiduciated  at  the  zero  of  the  open  circuit  voltage  profile. 
Subuact  15  ns  from  these  values  to  fiduciate  with  the  linear  rise  of  the  current  profile. 

Table  2 
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V.  IMPLOSION  STABILITY 


A.  Introduction 

Our  investigations  into  the  stabiUty  of  PRS  loads  have  shown  that  the  outer  radius 
is  limited  to  about  3  cm  if  the  load  is  a  thin  sheU.  Larger  radius  loads  are  Rayleigh-Taylor 
(R-T)  unstable  and  arc  prone  to  breakup  during  the  implosion.  Ihis  produces  a  poor  radiator 
at  stagnation.  We  have  found,  however,  that  the  stability  of  alternate  loads  mcluding  uniform 
fills,  multiple  shells,  and  tailored  density  loads  can  alleviate  this  constraint  on  radius. 

Z-pinch  implosions  can  be  characterized  by  three  major  stages:  pinch  formation  and 
initial  compression,  the  run-in  or  kinetic  energy  phase,  and  a  stagnation  or  radiation  phase. 
In  the  first  stage,  gas  breakdown  takes  place  and  an  imploding  shell  is  formed  as  the  material 
on  the  outside  edge  of  the  load  is  compressed.  Initially,  the  acceleration  is  slow  and  the 
internal  energy  produced  by  the  compression  is  much  larger  than  the  kinetic  energy  of  the 
plasma.  In  the  case  of  a  uniform  fill,  this  phase  and  the  next  phase,  the  run-in,  may  overlap 
but  there  is  still  a  period  where  these  criteria,  will  be  met  Due  to  the  low  acceleration, 
hydrodynamic  stability  is  not  an  issue  during  this  phase.  However,  any  initial  unstable 
activity  and/or  irregularity  in  the  initial  electrical  breakdown  of  the  material  may  seed  later 
instabilities.  The  duration  of  this  phase  varies  depending  on  the  driving  current  and  mass 
distribution  of  the  load.  During  the  run-in  ph^,  the  plasma  shell  picks  up  kinetic  energy  as 
it  accelerates  towards  the  axis.  If  the  initial  mass  distribution  of  the  load  is  limited  in  radial 
extent,  as  is  tire  case  for  foils,  multiple  wire,  and  high  Mach  number  gas  puffs,  the  mass  in 
the  shell  remains  fixed  during  the  run-in  and  the  compressed  shell  is  accelerated  uniformly. 
It  is  during  this  phase,  tiiat  the  Rayleigh-Taylor  instability  is  most  active.  This  is  compUcated 
in  the  case  of  a  uniform  fill  load  because  mass  is  continuously  being  accreted  onto  the  shell 
which  changes  the  acceleration  as  weU  as  the  thickness  of  the  shefii.  The  final  phase,  or 
stagnation,  occurs  when  the  accelerated  mass  reaches  the  central  axis.  At  this  point,  radial 
kinetic  energy  is  converted  to  internal  energy  of  which  a  portion  is  subsequently  radiated 
away  as  X-rays.  Large  distortions  caused  by  the  R-T  instability  will  decrease  tiie  X-ray  yield 
since  some  fraction  of  the  radial  kinetic  energy  will  be  transferred  to  axial  kinetic  energy,  and 
increased  internal  energy  as  well  as  enhanced  magnetic  field  penetration  will  limit  or  fully 
disrupt  the  compression  of  the  pinch.  In  this  report,  we  confine  ourselves  to  die  second 

phase,  or,  run-in  phase. 

The  Rayleigh-Taylor  instability  in  the  context  of  a  Z-pinch  is  essentially  the  same  as 
that  of  classical  hydrodynamics  except  that  the  support  of  the  heavy  fluid  is  by  the  magnetic 
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[iiTo] 


im 


raveW  by  tte  steU  durmg  Ite  implosion  and  the  amplification  of  an  inilial  pcitinbanon  is 
J  =  (lnA)2W4n  where  A  is  the  ratio  of  the  finsd  to  initial  ampUmdes.  If  an  amphficattoa 
famor  as  large  as  1.000  can  be  tolerated  and  a  wavelengd.  of  0.5  cm  is  chosen,  it  Mows  that 
te  distance  traveled  by  the  shell,  d.  is  less  than  2  cm.  Although  this  ts  an  ideahzed 
argument,  it  demonstrates  how  deleterious  this  instabiUty  is  to  the  integrity  of  large  radtut 


The  nimierical  model 
dimensional  (2-0)  magnetoh] 


which  we  use  to  simulate 


is  a  two 


D)  magnetohydrodyaamic  (MHD)  model  based  on  separate  ion  and  electron 


dEj  / dt+ (Pi  +  Q,)dV / dt  = -W  ■  qi  +  CV(Te -Tj) 


[4] 


To  fliese  equations  arc  added  Maxwell's  equations,  viz., 
VxB  =  47cJ/c 
V  X  E  =  -9B  /  cdt 


along  with  Ohm's  law. 


E  =  qJ-uxB/c  I'J 

In  Eqs.(l)  through  (7),  Pe4  are  the  material  pressures  for  electrons  and  ions,  Ee,!  are  the 
specific  energies,  p  is  the  density,  V  is  the  inverse  of  the  density,  arc  the  heat  fluxes,  C 

is  the  electron-ion  energy  exchange  term.  Prad  ^  radiative  power,  and  Tl  J  is  the  ohmic 
heating.  In  Eqs.  (2)  and  (4),  Qa  is  an  artificial  viscosity  used  for  numerical  stability  in 
regions  of  strong  compression.  Ihe  remaining  symbols  have  their  conventional  meaning. 

The  atr>Tnic  physics  model  contains  all  the  ground  states  and  an  extensive  manifold  of 
energetically  and  diagnostically  important  excited  states  distributed  throughout  the  various 
ionization  stages.  In  model  used  here,  the  radiation  transport  is  only  calculated  in  the  local 
zone.  Equation-Of-State  tables  are  constructed  from  this  model  and  used  in  the  dynamic 
calculations.  Table  lookups  for  Pfad.  Tc.  and  <Z>  based  on  internal  energy  and  density  arc 

employed.  The  tables  were  constructed  using  a  detailed  Collisional  Radiative  Equilibrium 
(CRE)  ionization  physics  model^.  The  simulation  uses  a  voltage  waveform  to  drive  a  circuit 
model  characteristic  of  the  generator  and  load  to  calculate  a  self-consistent  current 

During  the  implosion  phase  the  radiation  is  not  as  critical  as  it  is  during  the  final 
collsqise  phase.  Radiation  is  always  emitted  during  all  three  phases  of  the  implosion  but  it 
does  not  become  energeticaUy  dominant  until  the  plasma  is  reasonably  ionized  and  dense. 
However,  during  the  final  run-in  phase,  the  L-shell  radiation  from  moderate-  and  high-Z 
materials  becomes  energetically  important  and  can  account  for  a  substantial  fraction  of  the 
total  x-ray  yield.  The  circuit  model  is  equally  important  Our  model  uses  the  self-consistent 
plasma  resistivity  and  more  importantly,  the  plasma  inductance,  to  determine  the  magnitude 
of  the  driving  current  as  a  function  of  time.  The  initial  inductance  of  the  plasma  is  set  by  the 
radial  placement  of  the  return  current  rods.  For  these  simulations,  we  initialized  each 
calculation  with  a  load  inductance  of  2.5  nH.  The  lumped  inductance  and  resistance  of  the 
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Howth  rate  is  due  to  <te  fact  that  th<;  current  and.  consequently,  the  accelerauon  change 
continuously  during  the  implosion.  When  the  density  ratio  teaches  40  percent,  the 
rate  falls  helow  the  classical  growth  rate  values  detetmined  from  the  acceleration.  is 

roughly  me  same  point  at  which  flic  total  masses  of  the  initial  sheU  and  inside  fill  are  ^uaL 

Based  on  a  growth  rate  of  ZxlO*  sec-»,  the  time  to  amplify  an  initial  peituthanon  by  a  factor 
of  1  000  is  about  35  nsec.  In  these  simulations,  this  is  about  one-half  of  the  nm-in  time  an 
impiies  that  them  is  plenty  of  time  for  breakup  of  the  load  to  occur.  Growth  rates  which 

wetemeasutedatotherwavelengthscompaiedsiniilariywithclassicalgrowflirates. 

In  order  to  make  a  fair  comparison  between  these  two  kinds  of  load  designs,  we  have 

to  introduce  a  measuie  for  instahiUty  growth  which  takes  into  account  not  only  the  fastest 

growth  rate,  but  also  the  non-uniformity  of  perturbation.  This  is  ilhistrated  m  Figure  . 

where  the  fractional  mass  which  is  unstable  during  the  imploaon  is  plotted  as  a  function  of 

the  fill  percentage.  Here,  we  have  defined  the  fracdonal  mass  ratio  as  flie  mass  with  greater 

than  the  mean  aidal  kinetic  energy,  divided  by  the  total  mass  of  the  implosion.  The  pomts 
lepresent  time  averaged  values  measured  during  the  linear  growth  phase.  This  fraction  n, 
groatet  than  40  percent  for  an  imploding  thin  shell  For  a  uniform  fill  load,  however,  it « 
ro-t—H  to  below  25  percent.  Although  the  unifoim  fill  loads  experience  an  imtiaUy  higto 
acceleration,  the  instability  is  effective  in  disrupting  a  smaU  percentage  of  the  total  mass.  The 

implication  is  that  the  R-T  instability  is  mote  damaging  to  shell  loads. 

Implosions  from  7  cm.  5  cm,  3  cm,  and  2  cm  were  considered.  Hgnie  3  shows  the 
results  of  the  7  cm.  2-D  simulation  for  an  initial  gas  shell.  As  with  the  shnulaflons  shown  m 

Fig  1  the  initial  shell  thickness  was  0.8  cm.  The  wavelength  of  this  simulation,  and  m  the 

subsequent  simulations  discussed  below,  was  0.66  cm.  Density  contours  are  shown  at  0. 
66  82,  and  100  nsec.  A  smaU  peituibation  which  initializes  the  insubility  is  imposed  ear  y 
in  the  simularion.  The  instabflity  then  develops  untU,  at  100  nsec,  the  sheU  has  been 
disrupted  to  the  point  that  the  density  at  flie  center  of  the  shell  has  been  reduced  to  the 
background  vacuum  density.  At  this  point,  the  simuUtion  is  tenninated  doe  to  the  mal^ty 
of  the  MHD  fluid  code  to  handle  vacuum  over  die  entire  radial  span  of  the  gnd.  The  back 
edge  of  the  compressed  sheU  has  only  traveled  approxiniately  1.5  cm  before  this 
This  is  in  rough  agreement  with  the  simple  estimate  made  eariier  based  on  classical  R-T 
growth.  Because  the  shell  was  disrupted  before  significant  compressional  heating  at  on-axis 

stagnation  could  occur,  there  was  no  measurable  K-shell  yield. 

The  results  from  the  2  cm  radius.  2-D  simulation  arc  shown  in  Fig.  4.  In  the  second 
plot,  at  76  nsec,  the  mstability  has  pushed  material  away  from  the  center.  However,  the  load 
n^ss  reaches  the  axis  just  before  breaking  apart  at  84  nsec.  At  86  nsec,  a  spot  has  formed  at 
the  center.  Subsequently,  the  more  massive  plasma  at  z  =  0  cm  and  at  z  =  0.66  cm  reaches 
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the  level  of  axial 


al  5  predicts  that  for  this  case  the  time  inter™!  required  for  the  shock  wave  pertorbatioo  to 
change  sign  (after  that  the  perturbation  ampUtude  exhibits  decaying  oscdlanons,  remami^  a 
snaU  fraction  of  its  initial  value)  is  15  nsec.  This  is  very  close  to  the  osciUatron  tune  of 
nsec  measured  in  the  2-D  simulation.  This  osciUation  time  is  dependent  on  the  Atwood 
nnmber.  sound  speeds  bod.  upstream  and  downstream  of  the  impact  shock,  and  the  veloctty 

Similar  results  have  been  obtained  with  shell-on-sheU  implosions.  Quahtatively,  the 
same  conditions  exist  at  the  point  of  impact  for  both  shell-on-sheU  and  puff-on-shell 
implosions.  Differences  between  these  two  load  designs  are  due  to  the  acceleration  an 
stability  of  the  plasma  outside  the  interior  sheU  prior  to  the  time  of  collision.  Up  to  this  tune, 
an  outside  puff  will  behave  in  the  same  fashion  as  that  of  a  unifonn  fiU  implosion. 
However,  as  wi*  the  unifonn  fill  implosions,  the  penaity  for  the  enhanced  stability  of  e 
overan  configumlion  is  reduced  radial  kinetic  energy  of  which  a  direct  consequence  is  lower 

yield. 

D.  Density  Tailored  Loads 

Direa  drive  implosions  on  PRS  simulatois  with  long  risetimes  present  problems  for 
stabilixing  the  PRS  load  against  hydrodynamic  instabilities.  In  order  to  couple  sufficient 
kinetic  energy  into  a  load  which  is  not  too  massive,  the  inidal  radius  of  die  implosion  must 
be  larger  than  what  is  common  in  today's  PRS  experiments.  Since  the  amplification  of  the 
perturbation  is  independent  of  the  implosion  nine,  but  is  instead  dependent  on  the  distance 
over  which  the  runin  occurs,  a  mechanism  for  stability  must  be  found. 

If  the  load  is  tailored  properly,  the  acceleranon  of  the  snowplow  sheU  will  be 
reduced  during  the  initial  renin  phase.  Hie  convendonal  Rayleigh-Taylor  (R-T)  instability 
grewdi  rate  win  then  be  greally  reduced  or  even  stabilized.  He  kinetic  energy  itself,  on  the 
other  hand,  can  increase  during  this  phase  because  mass  is  being  accreted  onto  the  sheU 
while  the  velocity  remains  nearly  constant  That  is.  KE  =  l/2mv2  increases  beca^  the 
mass,  m.  increases  whUe  the  velocity,  v,  is  neatly  uniform.  Iliis  requires  that  the  density 
is  increasing  from  the  outside  to  the  inside  surface.  The  inside  radius  for  this  type 
of  load  is  set  so  mat  once  die  plasma  has  formed  a  thin  sheU  the  accelerared  renin  distance 
wffl  be  widiin  the  stabiUty  distance  of  roughly  2-3  cm  for  wavelengths  m  the  1/2  cm  range. 
Note  that  the  accelerated  renin  stable  renin  distance  may  be  greater  than  the  2  -3  cm  limit 
since  that  limit  is  based  on  a  zero  initial  velocity.  This  is  a  topic  which  wdl  be  mvesngated 

in  Ae  future. 
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growth  is  algebraic  or  exponential.  In  the  R-M  slow  growth  phase  this  quantity  should 
oscillate  around  a  fixed  value.  This  behavior  is  seen  in  Fig.  1 1  prior  to  150  nsec.  Thereafter, 
the  steep  rise  of  this  quantity  reflects  Ae  exponential  growth  during  the  acceleration  phase. 

E.  Conclusion 

At  large  diameters,  implosions  initiated  by  thin  gas  puff  shells  arc  extremely 
unstable.  The  distance  over  which  they  travel  before  being  disrupted  is  severely  limited  by 
the  large  R-T  growth  rate.  Based  on  the  measured  growth  rates  of  the  shell  implosions 
investigated  here,  they  are  limited  to  approximately  2-3  cm  run-in  distance.  This  effect  is 
reduced  for  uniform  fill  loads  where  the  imploding  shell  contains  only  a  fraction  of  the  total 

mass. 

In  terms  of  converting  kinetic  energy,  diin  sheUs  arc  die  most  efficient  load  designs. 
These  loads,  therefore,  place  an  upper  limit  on  radiative  performance  for  a  given  machine 
design  and  current  level  However,  as  we  have  shown,  these  loads  arc  most  susceptible  to 
instability  and  breakup.  Less  efficient  kinetic  energy  converters  are  uniform  fill  and 
structured  loads  such  as  puff  on  shell  sheU  on  ^eU,  and  tailored  density  profiles.  However, 
these  designs  are  mote  stable  and  can  be  placed  out  at  larger  initial  radii.  The  overall  result,  at 
least  for  the  conditions  studied  here,  is  that  these  types  of  loads  may  do  as  well  or  better  than 
thin  shell  loads  when  2-D  stability  is  taken  into  account  Uniform  fill  and  structured  load 
rfftsigns  also  have  the  advantage  of  easing  the  stiffness  of  the  machine  requirements  by  going 
to  large  radius  and  long  implosion  times. 
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Growth  Rate  (IxlO^  sec"l) 


Hgure  1 :  Growth  rate  for  the  Raylcigh-Taylor  instability  as  a  function  of  the  ratio  of  inside 
density  to  shell  density  for  a  wavelengdi  of  0.66  cm,  assuming  an  initial  shell 
radius  of  5  cm  and  sheU  thickness  of  0.8  cm.  The  shaded  region  represents  the 
range  of  classical  growth  rates  calculated  from  the  acceleration  measured  during 
the  run-in  phase. 
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function  of  die  ratio  of  inside  density  to  shell  density 


Z  (cm) 


Onsec 


Rgure  3:  DeK-ly  contouB  for  the  7  cm  thin  shcU  implosion:  (a)  0  nsec,  (b)  66  nsec,  (c)  82 
nsec,  and  (d)  100  nsec.  The  density  contours  are  spaced  linearly  and  normalized 

to  the  peak  value  at  each  time. 
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Radius 


'igure  3:  (Continues)  Density  contours  for  the  7  cm  thin  shell  implosion;  (a)  0  nsec,  (b)  66  mec,  (c) 
82  nsec,  and  (d),  100  nsec.  The  density  contours  are  spaced  linearly  and  normalized  to  the 

peak  value  at  each  time. 


Radius  (cm) 


(b) 

Figure  4:  Density  contours  for  flie  2  cm  thin  shell  implosion:  (a)  0  nsec,  (b)  76  nsec,  (c)  84 
nsec,  (d)  86  nsec,  and  (e)  96  nsec. 
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Figure  4:  (Continues)  Density  contours  for  the  2  cm  thin  shell  implosion:  (a)  0  nsec,  (b)  76  nsec, 
84  nsec,  (d)  86  nsec  and  (e)  96  nsec. 
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Radius  (cm) 
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Figure  5:  (Continues)  Density  contours  for  the  7  cm  uniform  fill  implosion;  (a)  74  nsec,  (b)  110  nsec^ 
(c)  126  nsec,  and  (d)  136  nsec. 
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Figure  8:  Density  contours  for  the  1  cm  puff-on-shell  implosion:  (a)  62  nsec,  (b)  68  nsec, 
(c)  74  nsec,  and  (d)  78  nsec. 
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Figure  8:  (Continues)  Density  contours  for  the  1  cm  puff-or 
(c)  74  nsec,  and  (d)  78  nsec. 
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Figure  1 


VI.  IMPLOSIONS,  EQUILIBRIA  AND  STABILITY  OF 
ROTATING  RADIATING  Z-PINCH  PLASMAS 


A.  Introduction 

The  idea  of  using  rotation  to  suppress  the  Rayleigh-Taylor  (RT)  instability  of  an 
imploding  plasma  goes  back  about  20  years  ago  to  the  LINUS  project  developed  at  the  Naval 
Research  Laboratory.  ^  Recently,  it  was  discussed  again  by  Rostoker  and  his  collaborators  ,  who 
proposed  to  use  spin  control  of  the  RT  instabiUty  in  high-density  Z-pinches.  The  technical 
feasibility  of  producing  a  spinning  plasma  column  which,  like  a  conventional  gas-puff  jet,  could 
be  used  as  a  Z-pinch  load  driven  by  a  pulsed  power  machine,  has  been  demonstrated  in  the 
experiment^.  Rotation  was  shown  in  this  experiment  to  produce  homogeneous  plasmas, 
smoothing  out  the  azimuthal  non-uniformities  like  those  typically  produced  in  gas-puff  Z- 
pinches  during  initial  gas  breakdown  and  observed  in  ref.  4. 

This  concept  of  controlling  the  RT  instability  is  aimed  at  stabilizing  the  final,  stagnation 
phase  of  implosion,  by  eliminating  the  acceleration  fliat  may  cause  instability.  Indeed,  the  inward 
motion  of  the  shell  coasting  to  the  axis  in  the  absence  of  any  external  forces  -  that  is,  with  its 
acceleration  identically  zero,  -could  be  stopped  and  reversed  simply  due  to  conservation  of 
angular  momentum,  see  example  in  Section  III.  (The  value  of  R  characteristic  of  radial 
deceleration  could  be  quite  large,  although  the  radial  component  of  the  acceleration,  R-Rd^ 
vanishes.)  Being  most  important  for  controlled  fusion  applications,  this  aspect  of  stability  does 
not  appear  to  be  a  major  issue  when  applications  like  plasma  radiation  sources  are  considered.  In 
this  context,  the  development  of  the  RT  instability  during  the  run-in  phase  is  one  mostly 
damaging  the  quality  of  implosions  (e.g.,  see  ref.  5).  The  dominant  mode  of  this  instability  is 
known  to  be  the  sausage,  m=  0  mode,  which  in  the  absence  of  an  external  axial  magnetic  field 
is  the  only  one  that  can  develop  without  bending  the  magnetic  field  lines.  The  filamentation  RT 
modes  with  nonzero  m  are  also  developing  in  imploding  plasmas,  but  their  growth  rates  are 
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smaUer  than  tee  of  the  sausage  mode  -  for  short  wavelengths,  by  factor  of  =0.643, 

see*’’.  Hence,  although  decreasing  the  amplitudes  of  the  initial  azimuthal  pertmbations,  either 
auongh  rotation  or  by  other  means  (like  increasing  the  number  of  return  current  rods)  would 
certainly  improve  the  quality  of  implosion,  it  is  unlikely  to  become  a  major  stabdizing  factor. 

Could  rotation  help  central  the  RT  instability  during  the  lun-in  phase?  This  is  one  of  the 
problems  addressed  in  the  present  paper.  The  earUer  results*  indicate  that  a  rotating  plasma 
column  is  mote  unstable  than  a  resting  one,  and  therefore  are  not  optimistic.  Nevertheless,  since 
the  presence  of  angular  momentum  is  shown  to  have  exactly  the  same  dynamic  effect  on  1-D 
motion  of  ideal  Z-pinch  plasma  as  the  presence  of  ftozen-m  axial  magnetic  flux,  it  is  natural  to 
inquire  if  it  has  also  a  similar  stabilizing  effect,  and  why  not,  if  it  does  not. 

Whether  or  not  the  Z-pinch  plasma  is  deliberately  spun  before  implosion,  its  deviation  of 
symmetry  is  likely  to  produce  some  angular  momentum  in  the  start-up  phase.  Being  conserved, 
and  increasingly  important  for  higher  compression,  non-zero  angular  momentum  can  strongly 

affect  possible  final  stages  of  the  implosion  In  particular,  spinning  can  change  the  flow  pattern 

drastically  if  in  die  absence  of  rotation  a  radiative  collapse  is  predicted  to  occur. 

This  paper  is  structured  as  follows.  In  SecSon  II,  the  Bennett  and  Pease-Braginskii’s 
steady-state  equilibria  (and  slowly  varying  “time-dependent  equilibria”)  of  rotating,  radiating 
cylindrically-symmetric  plasma  are  studied.  In  Section  m.  dynamics  of  rotating,  radiating  ideal 
plasmas  is  studied  analyticaUy.  In  Section  TV,  stability  of  a  rotating  Z-pinch  column  is  smdied 

using  a  very  simple  model.  In  Section  V,  we  conclude  with  a  discussion. 


B.  Bennett  and  Pease-Braginsldi  equilibria 

Let  us  start  with  the  equation  describing  radial  pressure  balance  in  a  cylindrically 


symmetric,  uniformly  rotating  2>-pinch: 

(l-bZ)^(nr)-l-  =  pnv 


(1) 
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wlere  Z  is  the  average  ion  charge,  12  is  (he  angular  velocity.  Multiplying  (1)  by  2®-"  and 

integrating  over  r  firom  0  to  we  obtain  the  Bennett  relatione 

f'2 

(1 + Z)  NT+  =  2^’ 

where  AT  is  the  ion  line  density,  T  is  the  average  temperature,  /is  the  moment  of  inertia  per  umt 
length  of  the  column.  This  could  also  be  presented  as 

(i+z)Nr='^(/*-/^).  ® 

where  is  a  characteristic  of  rotation  with  dimensionaUty  of  current: 

(kA)  =  3 1  (p ,  mg/cm)^  {V,  fcm/sec), 

V  being  some  average  velocity  of  rotation,  /z  =  A/  -  die  line  mass  of  the  column.  As  long  as 
die  current  I »  Iq,  Eq.  (4)  represents  the  conventional  Bennett  equilibnum.  For  instance,  with 
^=100  mg/cm  and  ^=100  km/sec,  we  find:  7^  =1  MA.  Thus  the  effect  of  rotation  could  be 
noticeable  if  the  characteristic  velocity  V  were  of  the  same  order  as  the  implosion  velocity, 
which  for  die  multi-MA  pulsed-power  machines  is  about  3  XlO’  cm/sec. 

Now  consider  the  simplest  equilibrium  profiles,  corresponding  to  flat  profiles  of 
temperature  and  current  density:  Tir)  =  T  =  const,  and  j{r)=  I /  TCR^  =const.  Then 
equation  (1)  is  easily  solved  for  n(r) .  Introducing  a  dimensionless  parameter 

(5) 

2(1 +z)r 

ciraracterisac  of  the  toteHou  (the  limite  X«1  aud  Z»1  coirespoud  to  low  aud  high 
importance  of  rotation,  respectively),  and  die  dimensionless  coordinate 

one  can  present  the  solution  as 
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(7) 


yt(yi)=  . — : — -  e5q)[-A<l  -1]^)]} 

^  jtR^  1-1  +  exp(-X)  ^ 

The  equilibrium  density  profiles  are  shown  in  Fig.  1  for  several  values  of  the  parameter  X. 
In  the  limit  A  0  the  conventional  parabolic  profile  n  =  (2N  1 7tR  XI  ”  reproduced. 

The  opposite  extreme  corresponds  to  uniform  density  n={N  where 

0(z)  is  the  Heaviside’s  step  function. 

Consider  now  an  infinitely  slow  variation  of  equilibrium  conditions.  The  current  I,  the 
pinch  radius  R,  temperature  T.  angular  velocity  and  the  parameter  X  may  vary,  whereas  the  ion 
line  density  N  remains  invariant.  Another  invariant  is  the  conserved  angular  momentum; 

M  =  JO  =  const. 


With  the  aid  of  (7),  we  can  express  the  moment  of  inertia  J  via  the  parameter  X,  and  substitute 


into  (8),  to  obtain 

[2-2A+^  -2exp(-Aff  _2^ - - 

^  ^  ■  X'[>l-H-exp(-A,)]  m,N 


(9) 


Now  consider  the  power  balance.  Following  Pease’  and  Braginskii",  let  us  assume  that 
Ohmic  heating  due  to  classical  Spitzer  resistivity  is  balanced  by  the  bremsstrahlung  radiative 
losses  from  the  optically  thin  plasma  column  (this  might  refer  to  a  Z-pmch  m  a  hydrogen  or 
deuterium).  The  power  density  of  Ohmic  heating  is  given  by  Qj  where  the 

Spitzer  expression  for  plasma  conductivity  has  been  used,  Ob  being  a  known 

constant  for  a  given  plasma  composition.  In  our  case  of  uniform  temperature  and  current  density, 
so  that  the  total  ohmic  power  per  unit  length  of  the  plasma  column  is 


Pi 


(10) 


The  power  density  of  bremsstrahlung  losses  is  given  by  Q-  =  where  Q,  is  also 

a  known  constant  Integrating  this  expression  with  the  aid  of  Eq.  (7),  we  find  the  total  power  of 
bremsstrahlimg  losses; 
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*  Af2X  -  3  +  4exp(- A)  -  exp(-2X)] 

Pr  =  2^  [Qrf^  -  27di^  [^- 1  +  e5q)(-/L)f 

Since  the  Ohmic  heating  is  supposed  to  be  balanced  by  radiative  coolmg, 

P  =  P 


(11) 


(12) 


after  some  algebra  we  can  transform  (12)  into 

I  2^  Z  (13) 

I -IpB^  3(21  -3+4  exp(-l)  -  exp(-2l)]  ‘ 

Here,  I  pg  is  the  conventional  value  of  the  Pease-Braginskii’s  current  in  the  absence  of  rotation 
(about  1.4  MA  for  deuterium  plasma),  and  the  left-hand  side  of  Eq.  (13)  tends  to  this  value  in  the 
limit  A  -+  0.  However,  Eq.  (13)  shows  that  with  rotation  taken  into  account,  the  power  balance 

(12)  is  possible  not  only  for  /  =  /«,,  but  for  any  current  greater  than  /«,,  provided  that  the 
value  of  X  is  sufficiently  high.  The  dependence  (13)  is  illustrated  by  Fig.  2. 

Suppose  we  maintain  an  almost  Pease-Braginskii  equilibrium,  with  I  very  close  to  Ipg , 
and  small  but  non  zero  A=  and  then  start  increasing  the  current  very  slowly,  so  tiiat  there  is 
enough  time  to  establish  the  power  balance  ( 1 2)  for  any  current  value  of/.  With  I  and  /  related  by 

(13) ,  we  can  express  R  via  X  from  (9).  This  relation 


R  (14) 

Ro-^Fa)^ 


where 

12-  2A+  ^  -  2exp(-A)f  ^,5^ 

^(2)  =  i+exp(-A)][2A-  3+  4exp(-A)  -exp(-2A)]  ’ 

together  with  (13)  expresses  in  a  parametric  form  the  relation  between  I  and  R  for  a  slowly 
varying,  “time-dependent  equihbiium  » 

Figure  3  presents  the  ratio  i?  /  i?o  versus  the  normalized  current  /  /  /«,  for  two  values  of 
the  parameter  0.1  and  0.01.  In  both  cases,  there  is  no  radiative  collapse:  the  radius  of  the  Z- 
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pinch  and  ftc  angnlai  frequency  of  to  rofaSon  adjust  in  snch  a  way  ftat  the  power  balance  is 
realized  with  finite  R  for  any  value  of/. 

Of  course,  higher  compression  is  needed  to  reach  the  balance  for  low  values  of  For 

instance,  with  Ai=0  0>  increasing  the  current  from  /  =  /»  to  /  =  3/„  should  be  accompanied 
by  an  almost  SO-fold  radial  compression,  that  is,  by  a  2.5Z10’.fotd  density  jump.  But  this  is  stffl 
smaller  than  a  10‘-fold  or  higher  density  increase  predicted  by  the  models  of  radiaHve  collapse, " 
involving  a  thousandfold  or  higher  radial  compression. 

Does  the  value  3o-0.01  conespond  to  an  unrealisticaUy  fast  rotation?  The  ratio  of  kinetic 
energy  of  rotation  £.  to  the  thermal  energy  of  the  Z-pinch  column  i«  a' 

.^<<1  equals 

^  -Xzil 
3 

which  means  only  2S10-'  even  for  g  =5/3  (and  for  radiating  Z-pinches,  lower  values  of  g  are 
appropriate).  A  small  deviation  from  ideal  cylindrical  symmetry  of  mass  and  current  distribution 
in  the  Z-pinch  load  would  be  quite  sufficient  to  produce  initiaUy  slow  )  rotation,  which 

would  nevertheless  be  sufficient  to  stop  the  radiative  collapse  long  before  the  ultra-high 
densities”'"'  are  produced.  We  will  show  that  the  same  conclusion  would  be  vald  if  dynamics 
of  compression  is  taken  into  account  more  accurately. 


C.  Self-similar  implosion  dynamics 

To  esdmate  the  degree  of  compression  which  could  be  achieved  when  an  optically  thin 
radiating  Z-pinch  column  rotates  tuound  its  axis,  we  use  exact  self-similar  solutions  of  ideal 

MHD  equations  with  radiative  losses: 
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(18) 


-w,  /y  1)QL 

(uV)r  +  (r  -  1)7V  •  u  =  -  (1+Z)n’ 

^-Vx(uxB)  =  0, 

(uV)u  =  -^V(nr) + X  B)  X  B , 

All  notation  here  is  standard.  The  power  density  of  radiative  losses  Q.  could  be  expressed  either 
by  the  above  fonnula  a  =  for  bremsstrahlung  radiative  losses,  or  by  Q.  = 

if  line  radiation  is  supposed  to  be  a  major  energy  sink.  We  did  not  include  the  Ohmic  heating  and 
magnetic  diffusion  terms  in  the  right-hand  sides  of  Eqs.  (18)  and  (19),  respectively,  because 
with  these  terms  we  would  not  be  able  to  separate  the  variables.  Physically,  this  means  that  we 
consider  currents  greater  than  flie  Pease-Braginskii’s  value  I pg,  so  that  radiative  energy  losses 

are  si^)posed  to  be  greater  flian  flie  Ohnuc  heating. 

Exact  self-similar  solutions  for  a  1-D  motion  with  uniform  (homogeneous)  deformation 
are  obtained  in  a  standard  way,  and  below  we  present  only  the  results.  (Details  and  references 
could  be  found  in^^’).  All  the  variables  are  presented  as  separable  function  of  the  self-similar 
coordinate  h  given  by  (6)  [vdiere  R  =  R(t)  is  the  time-dependent  radius  of  the  plasma  column, 
s  /J(f  =  0)]  and  of  dimensionless  degree  of  compression  a(t)  -  R(t)  /  Rc‘ 

Ur  =  Ro^'n>  u,=R^a-h\,  n  =  noN(7i)a-\ 

B,  =BeoHe(-n)oc-\  B,  =  B^H,(n)a-\  (21) 

where  i3o,  Ho.  %,  B,^,  B^  are  normalization  constants,  N(7]),  0(V)>  H,(Tl),aDd  H,(7]) 
are  self-similar  profiles  of  density,  temperature,  azimuthal,  and  axial  magnetic  field,  respectively, 
0  (f)  is  a  dimensionless  function  describing  loss  of  entropy  due  to  radiation. 

With  time  measured  in  units  of  =  (4;imino)^%  /  equation  of  motion  for  a(0 
could  be  presented  in  the  form 


87 


,  b_ 


(22) 


where  the  constants  are 


47til  +  Z)n,T^ 

Bio 


b=Qotl  +‘^* 


^00 


(23) 


The  parameter  p  has  the  conventional  meaning  of  ratio  of  thermal  pressure  (or  energy 
density)  to  the  magnetic  one.  The  parameter  p  includes  contributions  from  the  centnfugal  term 
and  from  the  axial  magnetic  field.  Since  the  two  positive  tenns  are  added,  the  effects  of  rotation 
and  of  axial  magnetic  field  on  the  radial  compression  are  quite  similar,  their  relative  importance 
being  determined  by  the  ratio  of  the  two  terms  in  the  expression  (23)  for  b.  To  explam  the 
similarity,  we  note  that  conservation  of  angular  momentum  implies =  const,  so  that  the 
centrifugal  acceleration  Q^RocR-\  Conservation  of  axial  magnetic  flux  implies 
B^R^  =const ,  so  that  external  force  per  unit  length  of  the  plasma  column  needed  to  contain  it 

is  ISk)ocR-\  Equation  of  motion  (22)  clearly  shows  that  the  presence  of  angular 

momentum  and/or  axial  magnetic  field  is  dynamically  equivalent  to  addition  of  some  ideal  gas 
with  g=2  to  the  Z-pinch  column. 


The  self-similar  profiles  iV(7|),  0(7} ),  (l?  ).  and  i/,(7]),  as  well  as  the  equation  for 

6  (t) ,  depend  on  the  mechanism  of  radiative  energy  losses.  For  bremsstrahlung  losses 


X  I-"* 

N{n)=[— 


©(??)  =  ~3 — . 


liJl 

3 


2 

J 


(24a) 


and 

where  dtaeneionless  parameter  £  =  (r  -  /  d  +  has  the 

meanmg  of  a  ratio  of  the  time  needed  to  radiate  away  all  the  thermal  energy  of  the  plasma 
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column  to  the  characteristic  time  U  (this  is  essentially  the  same  as  given  in").  For  line  radiative 


losses 


2V 


0(7] ) 


=1#. 


H.(»I)  =  ^(l8-24ij"  +  9r)T. 


(24b) 


and 

0=-e0-ia2(2r-3)^  (25b) 

with  e  =  (r  -  l)(4;i7n,)‘' /  (1  +  As  seen  from  Eqs.  (25a,  b),  the  self¬ 

similar  profiles  are  independent  from  the  parameters  b,  b,  and  e,  which  determine  the  time 
history  of  the  Z-pinch  motion.  In  Fig.  4a,  b,  the  dimensionless  profiles  of  A((7]),  0(7]), 
Hg  if]),  i/j(7]),  and  axial  current  density  7^(tj)  =  (1  /  !  di]\r\Hgir])\  =  rfNirj)  /  Hgirf) 

are  plotted  for  bremsstrahlung  and  line  radiative  losses,  respectively.  The  profiles  are  basically 
quite  similar,  with  current  density  more  sharply  peaked  near  the  axis  in  the  case  of  line  radiation. 

To  study  the  dynamics  of  the  Z-pinch  column,  Eq.  (22)  and  either  Eq.  (25  a)  or  (25b) 
should  be  integrated  with  the  initial  conditions  a(0)  =  1,  Ci:(0)  =0,  0(0)  =  1.  For  each  case 
the  solution  is  fully  determined  by  the  values  of  parameters  b,  b,  and  e.  To  start  with  the 
implosion  phase,  the  sum  ^  +^7  must  be  less  than  unity.  The  results  are  very  similar  for  both 
cases,  and  we  present  in  Figs.  5-6  only  those  calculated  for  line  radiative  losses,  y  =1.3,  e=0.03. 

Figure  5  is  plotted  for  the  case  of  no  rotation  and/or  axial  magnetic  field  (6=0,  p  =0.5). 
Then  the  solution  describes  a  radiative  collapse:  Z-pinch  radius  R  and  the  characteristic  value 
p/pl  proportional  to  the  exponential  of  its  entropy  vanish  after  a  number  of  bounces  in  the  same 
way  as  described  in".  With  nonzero  angular  momentum  and/or  axial  magnetic  field  (Fig.  6, 
plotted  for  6=0.1  and  P=0.4),  we  see  tiiat  no  radiative  collapse  does  take  place.  After  a  while,  all 
the  thermal  energy  gets  radiated  away,  just  as  in  the  case  of  Fig.  5.  Then,  in  contrast  with  Fig.  5, 
the  pressure  of  azimuthal  magnetic  field  is  balanced  by  the  centrifugal  force  and/or  by  the 


pressure  of  aziniuthal  magnetic  field.  With  no  inatability  that  conld  limit  (he  number  of  bounces 
taken  into  account,  the  final  state  of  the  Z-pinch  column  is  represented  by  undamped  oscillattons. 
Just  as  in  Section  fi,  no  radiative  coUapse  is  predicted,  unless  (he  both  the  angular  momentum 
and  the  magnetic  fimr  of  axial  magnetic  field  are  exactly  zero.  As  it  is  obviously  hard  to  produce 
real  Z-pinches  with  exactly  zero  angular  momentum,  we  conclude  once  again  that  the  radiative 
collapse  predicted  by  0-D,  1-D,  and  even  2-D  (r-z)  modeling  for  cylindrical  Zl-inches  with 
/  >  I„  might  be  a  consequence  of  unrealisticany  high  symmetry  imposed  on  these  simulations. 
Since  a  rotating  plasma  column  is  itself  unstable  (see*  and  Section  IV),  formation  of  compUcated 
3-D  patterns  appears  to  be  a  more  likely  result  of  its  evolution  than  an  essentially  one- 
dimensional,  cylindrically-symmetric  radiative  collapse. 


D.  Sansage  instability  of  a  uniformly  rotating  Z-pinch 


exhibit  kind  of  Rayleigh-Taylor  instability.  Let  us  find  out  whether  there  is  any  kind  of 
stabilization  due  to  rotation,  since  the  local  gravity  g(r)  =  Q^r  is  provided  by  the  centrifugal 


acceleration. 


For  our  illustrative  purposes,  we  use  the  ideal  MHD  model,  and  will  not  disciiss  its 
validity  here,  although  it  appears  to  be  applicable  to  the  case  of  solid  rotation  (see  ).  It  is 
convenient  to  use  the  non-inertial,  uniformly  rotating  frame  of  reference.  Then  in  the  unperturbed 
state  the  plasma  particles  are  resting,  and  their  motion  is  due  solely  to  perturbation.  The  time 
variable  is  thus  separated,  and  we  can  seek  solutions  proportional  to  exp(  J7  +  ikz) ,  where  G  is 
the  growth  rate,  k  is  the  perturbation  wavenumber  (the  exponential  factor  will  be  omitted  below). 
The  velocity  u  then  is  itself  a  first-order  perturbation,  which  could  be  conventionally  expressed 
via  the  displacement  of  a  plasma  particle:  u  =  d^  I  dt  =  r^.  Equations  (17)  and  (19),  which 
remain  valid  in  flie  rotating  frame,  are  integrated  in  a  conventional  way  to  give 


B<*>  =  Vx(|xB<°>), 


(27) 

(28) 


where  the  superscript  (1)  marks  the  first-order  perturbation.  In  the  equation  of  motion  (20)  we 
neglect  the  pressure  term,  and  second-order  term  in  velocity,  but  have  to  mclude  the  terms 
representing  the  Coriohs  force  and  the  centrifugal  force,  to  obtain 

P  2  p(0)|  =  ^[(V  X  B<‘>  )x  B<°>  +  (V  X  B<“>)x  B<*>]+  x  + 


We  find  from  (27)  that  for  m  =  0  perturbations  and  B^=0  in  the  imperturbed  state,  the 
only  nonzero  component  of  perturbed  magnetic  field  is  Using  (26),  (27),  and  expressing 
from  (28),  we  can  find  that 


(30) 
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where 


tht  dimtMiooless  eigenvalue  n  is  defined  as  V  =  2“  F  /  .  Subsdtotmg  (30)  into  the  r- 


component  of  Eq.  (29),  we  obtain  a  modified  Bessel  equation 


I  2"+  ^  J 

dr^  dr 


fora 


variable  7  defined  as  Y  =  -{k^r^  /  and  related  to  ;£^  by 


dY 


dr' 

Solution  of  Eq.  (31),  which  is  not  singular  at  r=0,  is  given  by 


I 

F(r)  = 


(v^-  +  8) 

V 


1/2 


(31) 


(32) 


(33) 


The  boundary  condition  on  the  perturbed  free  surface  of  the  pinch  states  continuity  of 
magnetic  pressure  there,  which  in  our  case  is  the  same  as  continuity  of  Be  (since  the  unperturbed 
velocity  of  rotation  is  parallel  to  the  magnetic  field,  no  electric  field  emerges  in  our  rotating 

frame  of  reference): 

Being  both  curl-free  and  divergence  free,  the  perturbed  magnetic  field  in  vacuum  could  be 
presented  as  V<P,  where  0  is  some  scalar  function.  Due  to  azimuthal  symmetry,  has  no 
azimuthal  component,  and  therefore  the  external  perturbation  of  B^  is  due  to  displacement  of 


the  boundary: 


dr 


(35) 


whereas  for  the  internal  pertuibation  we  find. 

(dB^  |,(ii)  =  Sf 


A=r 


rgo,  ^ 


^  R 


(36) 


Substituting  (35)  and  (36)  into  (34), 


and  using  Eqs.  (30),  (32),  (33),  we  find  the 


dispersion  relation,  which  is  a  transcendental  equation  for  F  : 


92 


(37) 


where 

and  fcR=-~(At^-8)'“,  (38) 

and,  once  we  are  interested  in  real-valued  solutions,  the  real  parameter  m  should  vary  between 
2V2  and  infinity. 

In  the  long-wavelength  limit  JtR  -4  0  [which  means  KT  2^,  see  (38)],  the  right- 
hand  side  of  Eq.  (37)  tends  to  a  finite  value  A-l-Jl,  so  that  the  real  root  K=Kq  =  3.1781  is 
found  as  a  numerical  solution  of  a  transcendental  equation 

2J2 

The  growth  rate  for  long  wavelengths  is  thus  proportional  to  the  wavenumber 

r  ^■^im  =  G.629nkR. 

Ko 

In  the  opposite,  short-wavelength  limit  /:i?  -4  o®,  the  leading  term  in  the  right-hand  side 
of  Eq.  (37)  is  12,  implying  that  k  should  be  of  the  same  order  (if  the  arguments  of  the 
modified  Bessel  functions  are  much  larger  than  their  orders,  then  the  ratio  of  Bessel  fimctions 
tends  to  unity).  Then  the  relation  I  2  is  equivalent  to 

r  ^Q(kRy'^=igky'\ 

where  ^  =X2^i?  is  the  local  gravity  acceleration  at  the  surface  of  the  plasma  column.  We  see 
that  the  classical  expression  for  RT  instabiUty  growth  rate  is  reproduced  here  in  the  short- 
wavelength  limit. 

To  calculate  the  dispersion  curve  F  (k)  shown  in  Fig.  7  (curve  1),  Eq.  (37)  was  solved 
numerically.  For  comparison,  the  RT  dependence  F  (k)  given  by  the  right-hand  side  of  Eq. 
(41),  is  plotted  there  as  curve  2.  No  stabilizing  effect  is  seen  for  short  wavelengtti.  As  for  long 
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wavelength,  the  linear  dependence  like  (40)  is  quite  typical  for  accelerated  Z-pinches  with 
distributed  current,  see  This  also  cannot  be  regarded  as  stabilization. 

Now  we  can  explain  why  rotation,  which  in  1-D  flow  produces  the  same  dynamic  effect 
as  the  axial  magnetic  field,  does  not  help  to  stabilize.  Indeed,  as  seen  from  Eq.  (22),  mtroduction 
of  both  of  these  terms  is  equivalent  to  addition  of  some  fluid  with  /  =  2,  effectively  shifting  the 
plasma’s  /toward  this  value.  But  variation  of  /within  finite  limits  is  known  to  produce  Uttle 
effect  on  sausage  instability,^^  with  the  dominant  RT  mode  in  the  short-wavelength  limit 
independent  from  /,  and  hence,  not  affected  at  all,  as  confirmed  once  again  by  Eq.  (41).  The 
stabilizing  effect  of  the  axial  magnetic  field  is  basically  two-dimensional.  This  stabilization  is 
due  to  the  energy  required  to  bend  the  force  lines  of  axial  magnetic  field,  which  is  represented  m 
the  dispersion  relations  by  the  tenns  proportional  to  There  is  nothing  to  match  it  m  the 

case  of  rotation. 


E.  Discussion 

We  have  shown  that  plasma  rotation  is  not  likely  to  help  stabilize  the  most  dangerous 
sausage  modes  of  the  RT  instability  of  imploding  Z-pinches  in  the  run-m  phase.  Although  it 
might  be  helpful  in  stabilizing  filamentation  modes  (m>2),  (see  Ref  21)  as  well  as  m 
decreasing  the  initial  amplitudes  of  the  corresponding  perturbations,  since  these  are  not  dommant 
instabiUty  modes,  it  is  not  likely  to  improve  overall  implosion  stability.  Consequently,  the  idea  of 
making  a  rotating  Z-pinch  load  shows  Uttle  promise  for  applications  like  plasma  radiation 
sources,  unless  a  major  stabilizing  effect  of  sheared  rotation  (not  studied  here)  on  the  sausage 
instabiUty  modes  is  discovered. 

In  real  imploding  Z-pinches,  however,  some  angular  momentum  is  always  present,  even 
if  no  efforts  are  made  to  produce  it  deUberately.  We  have  demonstrated  above  that  even  small 
initial  angular  momentum  is  sufficient  to  stop  radiative  collapse  for  a  current  above  the  Pease- 
Braginsldi  value  long  before  the  ultra-high  densities  are  achieved  (for  instance,  starting  rotational 
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kinetic  energy  of  about  10*^  of  the  fliennal  energy  is  enough  to  limit  radial  compression  of  a 
quasi-steady  plasma  column  carrying  current  /  =  3/„,  at  the  level  below  50).  This  is  why  die 
radiative  collapse  of  a  cylindrical  plasma  column  predicted  for  I  >  Ipg  might  be  a  consequence 
of  unrealistically  high  symmetry  imposed  on  analytic  studies  and  numerical  simulations, 

including  even  2-D  ones. 

Both  trapped  flux  of  axial  magnetic  field  and  nonzero  angular  momentum  are  capable  of 
stopping  radiative  collapse.  Although  their  dynamic  effects  in  ideal  MHD  model  are  shown  to  be 
similar,  the  mechanisms  of  their  dissipative  losses  are  quite  different.  The  axial  magnetic  field 
can  simply  diffuse  through  plasma  into  vacuum,  and  it  could  easily  be  dissipated  if  the  plasma 
radius  is  sufficiently  small.  In  contrast  with  this,  angular  momentum  could  only  be  lost  m  axial 
direction,  either  through  ejection  of  mass  carrying  momentum,  or  through  viscosity.  The  former 
mechanism  is  evidently  more  effective.  Hence,  the  actual  high  compression  of  radiating  plasmas 
in  Z-pinches,  which  leads  to  formation  of  bright  spots,  being  necessarily  accompanied  by  loss  of 
angular  momentum,  should  produce  not  cylinder-like  regions  of  dense  plasma  (U-shaped  necks), 
but  rather  point-like  ones  (V-shaped  necks).  This  seems  to  be  in  agreement  both  with  the 
available  theoretical  data,  and  with  earlier  theoretical  results^  based  on  different  assumptions. 
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R.  p/p^ 


t 

Fig.  6.  Same  as  in  Fig.  5  in  the  presence  of  rotation 
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r/a 


kR 

Fig.  7.  Growth  rate  Fexpressed  in  units  of  angular  velocity  Q  vereus  Ai?:  1,  numerical 
solution  of  Eq.  (39);  2  -  classical  RT  growth  rate  (41). 
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x-ray  emission  spectra  from  the  K-sheU  provide  valuable  information  wtth  wtah  to  asse^ 
plasma  conditions  and  X-ray  emission  performance  in  Z-pinch  plasmas.  Th^  soft  “ 
„«tinely  used  to  diagnose  K-ions-  for  both  laboratory  and  astrophys.^  pla^as.  Smce  r^o 

involving  Hand  He  lines  are  so  useful  in  obtaining  diagnostics  mlow-Zpmches,  ratio 

for  F  and  Ne-lto  systems  should  provide  useful  X-ray  diagnostics  for  higher  atom.c  number!  ) 
Z-pinch  plasmas.  Since  the  L-sheU  ions  are  many  electron  systems,  then  X-ray  " 

in  structure  due  to  a  large  number  of  msonance  and  sateUite  lines,  many  “  .Tm  LTta 
one  another.  Identiflcation  of  these  lines  for  diagnostic  calculadons  pose  chaUenges  that  cm 
addressed  and  hence  diagnostics  r«earch  involving  I.sheD  ions,  while  not  as  ~  ^ 

involving  K-shell  ions,  offers  potency  more  rewards.  Efforts  have  already  been 
the  I^shen  diagnostics  of  Ne-like-  ions  due  to  its  closed  shell  conaguranon  for  the  ground  state. 

The  intensity  rahos  of  Ne-  and  F-like  emission  lines  are  sensidve  to  electron  temperatme 

density,  and  rates  of  ionirahon  and  different  line  rahos  need  to  be  used  to  infer  “ 

_ _ _ _  TU.  rahos  of  two  allowed  lines  of  the  same  ton  generally  cannot  be 


satellite  lines.  The  intensity  ratios  of  two  resonance  lines 
■t:  H  .rnod  source  for  inferring  the  state  of  ionization  stages 


ion  suges  can  generate  a  temperature-dependent  expression  for  the  relative  ion  rracuoi^  xu. 
"nTs  of  the  resonance  lines  of  two  neighboring  ions  is  aHo  density  se^uve  We  ^ 
investigate  the  various  line  ratios  as  diagnostic  tools  in  both  thne-independent,  colltsional-radiattv 

equilibrium  (CRE)  and  transient  plasmas. 

In  Z-pinch  plasmas  the  principal  excited  levels  are  populated  by  electron  impact  exdration 
from  the  ground  state,  and  depopulated  by  radiative  decay.  When  these  processes  are  much  fasmr 
!  Changes  in  ie  plasma,  steady  state  equations  can  adequately  de^nbe 


a.e  plasmas.  For  mid-Z  plasmas  that  cao  be  made  to  radiate  in  the  K-sheU  however,  iomzation 
and  recombination  relaxation  times  in  the  I^sheU  may  be  longer  than  the  plasma  hydrodynamic 
times  and  time-dependent  or  transient  equations  must  be  used  for  ionixadon  calculadons.  But,  m 
order  to  inter  plasma  characteristics  using  time  dependent  calculations  and  then  compare  them  to 
experimental  observations  one  will  require  time  resolved  spectra  that  are  some  times  not  available. 
Satellites  of  multiply  charged  ion  resonance  lines  in  time  integrated  spectra  might  however  be 
used  as  reliable  ionization  gauges  provided  these  line  ratios  can  be  calculated  usmg  detailed 
and  accurate  time  dependent  calculations.  If  fomrc  experimental  endeavors  are  geared  towards 
generating  spectrally  infonnative,  spatially  and  time  resolved  data,  then  it  will  need  to  be  analyzed 
in  accordance  with  accurate  theoretical  predictions. 

This  present  diagnostic  work  involves  the  strong  resonance  and  satellite  lines  of  transitions 
from  the  n=3  levels  of  F-  and  Ne-like  selenium.  The  n=3  F-like  and  Ne-like  states  of  selenium 
are  non-statistically  distributed  in  the  density  region  of  diagnostic  interest,  10“  <  <  10  ,  and 

this  provides  the  basis  for  plasma  diagnostics  using  line  emissions  from  these  states.  Figures  1(a) 
and  (b)  show  the  I^shell  spectrum  of  selenium  and  krypton  under  low  spectral  resolution  so  that 
the  sateUite  structures  are  not  resolved.  Since  the  Ne-like  ground  state  has  a  closed  basic  shell 
structure  and  the  ground  state  configuration  2a’2p«  has  only  one  level  ^5o,  the  short-wavelength 
spectrum  is  simple  having  only  seven  allowed  resonance  transitions  and  these  lines  along  with  their 
satellites  can  be  reliably  identified  and  detected  provided  experimental  observations  are  earned  out 
with  high  spectral  resolutions.  A  high  resolution  X-ray  spectrum  for  a  mid-Z  Ne-like  ion  is  shown 

in  Fig.  1(c). 

For  F-lflre  ions,  on  the  other  hand,  the  spectrum  is  quite  compUcated.  Here  are  34  allowed 
transitions  of  the  2p'.2p*3d  type  alone  between  the  2p.3d  and  2.-3p  Ne-like  iines.  For  reliabie 
diagnostic  work  using  individual  lines,  it  is  important  that  these  resonance  hnes  be  isolated  and 
properly  identified  as  individual  lines.  The  satellites  of  these  F-like  2p-3d  and  2p-3a  resonance  lines 

haveavery  large  number  oflines  whose  individual  detection  with  absolute  certainty  wOl  be  difficult 

even  with  extreme  high  experimental  spectral  resolution.  Thus  for  the  F-like  re«.nance  hues  we 
the  sum  of  overlapping  lines  originating  from  the  same  configuration  as  one  smgle  fine 
cluster  For  example,  all  the  lines  for  the  2p-3a  resonance  transitions  are  summed  up  and  is  treated 
as  a  single  transition  and  the  same  treatment  is  Mowed  for  2p-3d  and  2.-3p  transitions.  We  do 
flie  same  lumping  when  we  consider  the  Ne-like  satellites  to  these  F-like  resonance  lines.  These 
20-34  and  2p-3«  line  dusters  for  F-like  and  Ne-like  selenium  are  non-overlappmg  and  distinct  and 
therefore  the  intensity  ratios  involving  these  lines  as  references  are  reliable  for  plasma  diagnostics. 
We  analyze  the  individual  fine-structure  lines  of  the  Ne-like  spectra  inorder  to  obtain  an  accurate 
sum  When  we  calculate  the  intensity  ratios  of  F-like  to  Ne-like  resonance  lines,  we  do  the  same 


lumping  for  the  Ne-like  lines  by  summing  those  individual  ones. 


B. 

We  have  developed  a  prototype  atomic  model  specifically  for  the  fiuorme  and  neon  ionization 
stages  that  is  necessary  for  analyzing  highly  resolved  X-ray  line  emission  from  these  tomzanon 
Mages.  This  model  can  be  nsed  in  the  analysis  of  wire  implosions  tf  selemum  is  coated  on  as 
a  diagnostic  element  or  it  can  he  nsed  in  the  analysis  of  gas  pnff  implosion 
as  an  impurity  element  file  energy  level  diagram  of  our  model  is  shown  m  Fig.  •  ®  ' 

like  model  consisted  of  31  states  including  the  ground  level,  26  n=3  fine-structure  of 
2p'3a  2p‘3p  and  2p=3d  configurations.  3  lumped  inner-sheU  excited  configuranons  of  2s2p  3s, 
2s2p'3p  and  2s2p‘3d  and  a  lumped  »=4  level.“  We  have  updated  our  F-like  tomzation  stage  by 
including  the  fine-stracmm  levels  for  the  «=3  singly  excited  states.  Thus  the  F-bke  stage  conl^  a 
lumped  ground  state,  the  A»=0  excited  state.  57  fine-structure  levels  of  the  2p  3s,  2p  3p,  ^d  2p  3 
configurations,  3  lumped  iuner-shen  excited  configurMions  of  2»2p=3s,  2s2p‘3p.  and  2s2p  3d  ^d 
two  lumped  «=4  states  For  the  Na-like  n=3  excited  states,  all  fine  strucmte  levels  of  the  2p  3 
as  weU  as  doubly  excited  2p‘3i3/'  and  2s2p*3I3l'  states  are  included.  We  have  also  considere 
fine-structure  levels  of  the  double  excited  2p‘3/3i'  and  2s2p‘3l3i'  levels  of  the  Ne-like  “ 

ti,e  calculation  of  satellite  lines  formation.  The  fluorine-like  and  neon-like  models  are  embedded 
in  a  detaUed  configuration  accounting  (DCA)  atomic  model  of  selemum  from  which  the  tomzation 
balance  is  calculated.'’  The  n=3  singly  excited  states  for  both  F-  and  Ne-like  ions  are  fiiUy  coupled 
to  the  ground  states  and  among  themselves. 

The  Mowing  linefonnation  mechanisms  are  considered  in  our  calculations.  For  the  resonance 

lines  the  most  dominantcontribntions  come  from  direct  electron-impact  exctanonfrom  the  ground 

state  We  have  included  direct  collisional  excitations,  deexcitations  and  radiative  decay  from  or 
to  the  ground  states,  among  the  n=3  excited  states  and  cascades  from  the  n-4  state  of  Se  an 
Se’«  In  this  calculation  we  have  aUo  included  resonance  excitation  to  the  n-3  excited  states 
of  Ne-like  selenium.  Collisional  ionization  from  the  n=3  states  of  both  F-  and  Ne-hke  tomzation 
stages  are  included  in  this  work.  The  excited  states  are  substantially  populated  by  DR  an  we 
have  included  a  detailed  state  specific  contribution  of  DR  for  recombitmtion  from  the  ground  and 
A„=0  state  of  Se-  and  Se25.  to  the  ,=3  singly  excited  states  of  Se-  and  Se-  respectively 
The  contributions  from  the  DR  sateUites  due  to  high  Rydberg  levels  that  are  unresolved  from 

resonance  lines  have  important  consequences''  and  they  must  he  ufcn  into  -omit  -n  -dm 
,0  compam  the  hue  ratios  with  experiments.  Previous  diagnostic  models  were  unable  to  predmt 
accurately  because  they  lacked  these  necessary  theoretical  considerations  m  their  calculations  and 

experiniciital  refiiiciiients. 
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Satellite  lines  are  formed  on  the  long  wavelength  sides  of  the  resonance  lines  if  the  spectator 
electron  is  in  the  n=3  states  of  the  doubly  excited  states  that  are  formed  due  to  the  DR  process. 
As  mentioned  before,  for  the  n  >  4  states,  the  satellite  lines  formed  are  often  blended  with  the 
resonance  lines.  For  the  formation  of  the  Ne-  and  Na-like  satelUte  lines  we  have  included  the 
results  of  detailed  calculations  of  DR  from  the  ground  states  of  Se“+  and  Se^*+  states  to  the  3/31' 
doubly  excited  states  of  each  ion.  TTiese  satellite  lines  can  also  be  formed  due  to  an  alternative 
mechanism  of  direct  inner-shell  excitation  (ISE)  of  an  electron  of  the  ground  or  the  An=0  state  of 
the  ion.  We  have  included  the  contribution  of  ISE  for  the  formation  of  Na-like  satelUte  lines.  The 
sodium-like  and  neon-like  ground  states  are  independent  source  terms  and  thus  the  populations  of 
the  upper  levels  for  the  satelUte  lines  due  to  DR  from  the  ground  state  and  ISE  from  the  Na-like 
states  are  added  as  separate  contributions. 

The  atomic  data  for  this  calculation  was  either  calculated  or  obtained  from  very  detailed  and 
sophisticated  atomic  codes.  The  atomic  structure  data  such  as  the  energies,  osciUator  strengths 
and  the  radiative  rates  for  aU  the  64  F-Uke  levels  were  calculated  from  the  atomic  code  developed 
at  Penn  State  University  by  H.  L.  Zhang  and  D.  H.  Sampson.  Electron-impact  colUsion  strengths 
for  excitation  from  the  ground  to  aU  the  excited  states  and  the  colUsional  coupling  among  aU  the 
excited  F-Uke  levels  are  calculated  by  using  the  relativistic  distorted  wave  code  (DFW)  of  Sampson 
et  al}^  and  the  electron  impact  ionization  cross  sections  for  ionization  from  the  ground  as  well 
as  as  excited  F-Uke  states  were  also  calculated  from  the  ionization  DFW  codes.  DR  rates  for 
recombination  from  O-  to  F-Uke  states  were  calculated  by  us“  using  Cowan’s  code.*°  For  Ne-Uke 
selenium  coUision  strengths  for  excitation  from  the  ground  state  were  obtained  from  Zhang  et 
and  the  colUsional  couplings  among  the  n=3  levels  were  obtained  from  Hagelstein  and  Jung.”  The 
colUsional  ionization  cross  sections  from  the  n=3  Ne-Uke  levels  were  obtained  from  Golden  er «/.” 
and  Moores  et  al.^*  The  DR  data  for  recombination  from  F-  to  Ne-Uke  levels  were  calculated  by 
us.“  TUe  excitation  and  ionization  rate  coefficients  were  obtained  by  numencaUy  integrating  the 
colUsion  strengths  or  cross  sections  over  a  MaxwelUan  distribution. 

The  power  emitted  per  unit  volume  for  the  F-Uke  or  Ne-Uke  resonance  lines  are  given  by: 

i 

and  the  power  emitted  per  unit  volume  in  a  Ne-Uke  or  Na-like  satelUte  line  is: 

Ar 

Ps  =  NelNp/NeCd  +  Nfre/Na^e]  ^  ^ 

Where  a,  is  the  additional  contribution  due  to  unresolved  satelUte  lines  emission  and  the  power 
emitted  for  the  Ne-Uke  satelUte  lines  do  not  include  the  contributions  due  to  ISE  in  this  present 
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calculation.  «  and  N;  are  densities  of  the  ground  and  eacited  sates  of  *e 

electron  density,  Cb  and  Cdb  are  the  exciation  and  DR  rate  coeffiaena,  0  -  •) « 

decay  aB  from  level  f  to  level  i  and  is  the  energy  for  this  transition,  the  wulatton  dens^ 

of  the  exciad  sates  aat  are  analyxed  for  both  CRE  and  time  dependent  calculahons  and  the  tm^e 

histories  of  electron  and  ion  amperatures  are  imported  into  our  analysis  from  the  coupled  rate 

equation  calculations. 


C.  Resolte 

In  this  present  diagnostic  calculations  we  considered  the  2p-3s,  2p-3d.  and  2s-3p  reronance 

cl„sarsandd.eirsaamtesforbothF-lilaandNe-likeselenitnn.Theiine-strnctu^ 

for  Ne-like  selenium  are  identihed  in  Table  1.  To  test  the  ionization  calculatton,  a  prototyptcal 
case  of  rapid  heating  and  cooling  was  investigated.  The  time  histories  of  the  electron  and  ton 
amperatures  are  shown  in  Fig.  3.  The  rapid  time  variation  of  the  electron  amperamres  wtU 
affect  the  ionization  stages  of  selenium  and  alar  the  strengths  of  recombination  vs  exettation  ^d 
ionization.  In  fact,  when  the  plasmas  bum  through  the  L-shell  to  the  K-shell  quia  mp.  y, 
wUl  have  significant  effects  on  the  line  ratios.  Figure  4  shows  the  tune  vananon  of  the  three  - 
^Ne-like  Inance  line  ratios  for  the  electron  density  of  5.95  x  1(FW».  The  am^ramre 
sensitivity  of  the  F  a  Ne  resonance  line  ratios  for  two  different  denmties  are  shown  m  Ftgs.  5 
and  6  As  can  be  seen  from  these  two  graphs,  the  inansity  ratio  of  the  2p-3.  resonance  hne  ts 
insensitive  to  tempemture  variations  for  eBctron  amperatures  r.>400  eV  and  thus  ts  not  a  good 
aol  for  amperamre  diagnostic.  The  density  de,andences  of  da  F  a  Ne  r«onance  h^  mtios  ^ 
be  seen  in  Figs.  7  and  8.  The  density  sensitivities  of  ratios  of  resonance  toes  for  both  F  -hke 
and  Ne-like  selenium  is  shown  in  Fig.  9.  In  the  experimental  observations  of  X-y  r« 
lines  from  Ne-like  bromine.’  the  line  ratios  were  found  a  be  density  sensttive.  At  low  denstti« 
ae  majority  of  the  exciad  populations  are  in  the  2p’3a  and  2p’3p  staas  and  most  of  the  manstty 
rre  2P-3.  resonance  lines  cmne  ffom  cascade  a  the  3a  saas  from  ae  higher  «c.a  sate. 
As  the  density  increases,  some  of  the  higher  exciad  metastable  staas  depopulaa  by  cotaonM 
deexcitation  a  lower  lying  exciad  staas.  This  causes  the  2p-3./2.-3p  resotmee 
density  dependent  This  density  dependency  due  a  depopulation  of  metastab le  exaad  staa 
tire  line  ratios  of  2p-3d/2»-3p  and  2p-3a/2a-3p  for  both  F-like  and  Ne-hke  selenium  can  seen 

from  Fig.  9. 

The  line  ratios  of  saallia  a  the  resonance  lines  for  both  F-lito  clusar  Itos 
„  lines  provide  extremely  good  amperamre  diagnoses  Ft^  ^ 

line  ratios  for  F-like  and  Ne-like  selenium  for  an  electron  density  of  10  cm  .  ^  can  be  seen 

tomFig  IktheconUibutionatheinansityofthesaallialinesduetoISEissigmacantforsome 
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sttong  lines  compared  to  some  others.  The  effect  of  ISE  on  the  formation  of  the  sateUim  hnes  a^ 
hence  on  the  line  ratios  become  more  important  with  the  increase  of  plasma  heatmg.  Hus  can 
seen  from  Fig.  12  where  we  show  these  line  ratios  for  a  transient  calculation.  The  ratio  of  the 
sateUite  lines  whose  upper  leveU  are  formed  including  both  DR  and  TSE,  to  the  resonant  hnes 
for  Ne-lihe  selenium  depend  strongly  on  the  ratio  of  Ne-like  to  Na-like  ion  populanons.  Thus  as 
the  plasma  bums  through  the  ionization  stages,  the  depletion  of  Na-like  selenium  tons  compar^ 
to  that  of  the  Ne-like  ions  win  have  important  consequences  on  the  hne  tanos  and  hence  on  e 

diagnostics  of  plasma  parameters.  The  ratio  of  the  sum  of  the  satellites(»>3)and  that  of  the  sum  of 

the  same  sateUite  totesonancelinescan  be  used  toobtainafractionalpopulation  abundance  rano  of 

two  ionization  stages  in  selenium.  ’Ihe  maximum  plasma  electron  temperature  can  be  determme 
from  ±e  dielectric  component  of  a  sateUite  to  the  excitation  component  of  a  resonance  line  for 
both  F-like  and  Ne-like  selenium  in  time  dependent  calculations.  Hie  time  variation  of  the  ratios 
of  sateUite  to  resonance  lines  for  Ne-like  selenium  is  much  more  sigraficant  than  those  for  F-  e 
selenium  For  the  F-Uke  calculations  we  consider  the  clusters  of  Unes  as  a  single  line  and  we  so 
do  not  include  any  contributions  of  ISE  of  the  Ne-like  sutes  to  the  intensity  calculations  for  the 
sateUite  lines  and  the  for  some  cases  the  contributions  of  ISE  increase  the  line  ratios  by  orders  of 

magnitude. 


D.  Summary 

We  have  considered  a  Maxwellian  electron  distribution  for  this  present  calculations.  Howew. 
it  is  possible  that  intense  heating  wiU  generate  non-MaxweUian  distributions  which  are  predicted  to  ■ 
have  non-negUgible  effects  on  the  line  ratios  and  hence  on  the  diagnostics  of  these  strongly  heated 
plasmas.ti“.“  To  complete  these  calculations,  we  wUl  include  the  effects  of  non-Maxwelhan 

distributions  in  the  future. 


It  is  apparent  that  in  future  high-current  driven  Z-pinch  experiments  involving  medium  to  high 
Z  plasmas,  rapid  ‘1.0™  through”  the  L-sheU  is  required.  Thus  an  accurate  and  detaded  analysts  of 
tite  line  emission  in  the  L-sheU  spectra  of  these  plasmas  wfll  be  extremely  relevant  ^  a 

nnderstandingof  the dynamicsof  medium  to highZmaterialpinchesthatarecom^d of  dtferent 

chargedstates.HiesediagnosticsmdiesfordensitiesachievedweUabovelO  , where excit^sta 

are  coUisionaUy  coupled,  it  is  essential  to  have  an  atomic  model,  such  as  oms  that  mcludes  al 
the  dominant  atomic  processes  to  aU  the  important  excited  levels  and  a  detailed  leve  ac^unting 
cannot  be  simpUfied  by  an  averaging  process  such  as  used  in  the  average  atom  model,  wbch 

inhibitauthehnpormntandrelevantatomicdynamics.  litis  point  is  emphasized  by  Ftg.  13wbch 

aiusnates  graphicaUy  how  stmctutaUy  different  the  n=3  states  of  F-  and  Ne-hke  tom  are. 
present  study  will  also  help  us  to  predict  how  long  the  plasma  remains  m  the  l^sheU  and  how 
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fast  wffl  the  plasma  ionize  and  bum  through  the  K-sheU  states.  We  can  conclude  that  the  analysis 
presented  here  can  predict  the  energy  absorbed  per  unit  mass  of  the  plasma  and  the  amount  of  mass 
participation  in  the  L-sheU  plasma.  Such  mass  inferences  provide  a  good  indication  of  the  quality 
of  implosion  and  of  the  X-ray  production  in  Z-pinch  experiments.^ 

References 

1.  J.  P.  Matte,  J.  C.  Keifer,  S.  Ethier,  M.  Chaker,  and  O.  Peyrusse,  Phys.  Rev.  Lett.  72, 1208  (1994). 

2.  M.  C.  Coulter,  K.  G.  Whitney,  and  I.  W,  IhomhUl,  J.  Quan.  Spectre.  Radiat  Transf.  44, 443 
(1990). 

3.  V.  A.  Boiko,  S.  A.  Pikuz,  and  A.  Ya.  Faenov,  J.  Phys.  B  12, 1889  (1979). 

4.  A.  Zigler  et  al.,  Phys.  Rev.  A  45, 1569  (1992). 

5.  W.  H.  Goldstein,  R.  S.  Walling,  J.  Bailey,  M.  H.  Chen,  R.  Fortner,  M.  Klapisch,  T.  PhilUps,  and 
R.  E.  Stewart,  Phys.  Rev.  Lett.  58, 2300  (1987). 

6.  S.  Ya.  Khakhalin,  B.  A  Bryunetkin,  I.  YU.  Skobelev,  A.  Ya.  Faenov,  J.  Nilsen,  A.  L.  Osterheld, 
and  S.  A.  Pikuz.  JETP  78, 633  (1994). 

7.  J.  Bailey,  R.  E.  Stewart,  J.  D.  Kilkenny,  R.  S.  Walling,  T.  Phillips,  R.  J.  Fprtner,  and  R.  W.  Ixe, 
J.  Phys.  B  19, 2639  (1986). 

8.  D.  A.  Liedahl,  S.  M.  Khan,  A.  L.  Osterheld,  W.  H.  Goldstein,  in  High  Resolution  X-Ray 
Spectwsepy  of  Cosmic  Plasmas,  edited  by  P.  Georenstein  and  M.  Zombeck  (Cmabndge: 
Cambridge  Univ.  Press  1990)  lAU  Colloq.  115,  p.  49. 

9.  A.  H.  Gabriel,  Mon.  Not.  R.  Aston.  Soc.  160, ,99  (1972). 

10.  A,  H.  Gabriel  and  T.  M.  Paget,  J.  Phys.  B  5, 673  (1972). 

1 1.  A.  Dasgupta  and  K.  G.  Whitney,  Phys.  Rev.  A  42, 2640  (1990). 

12.  K.  G.  Whitney  and  M.  C.  Coulter,  IEEE  Trans,  on  Plasma  Sci.  16, 552  (1988). 

13.  J.  C.  Raymond  and  B.  W.  Smith,  Astrophy.  J.  306, 762  (1986). 

14.  E.  V.  Aglitski,  V.  A  Bioko,  O.  N.  Krokhin  et  al.,  Kvant  Elektron.  1, 2067  (1974). 

15.  V.  A.  Boiko,  A  Ya.  Faenov,  and  S.  A.  Pikuz,  J.  (JuanL  Spectr.  Radiat  Transf,  19, 11  (1978) 

16.  P.  G.  Burkhalter  and  D.  J.  Nagel,  Phys.  Rev.  A 11, 782  (1985). 

17.  J.  F.  Seely,  T.  W.  Phillips,  R.  S  Walling  et  al.,  Phys.  Rev.  A  34, 2942  (1986). 

18.  D.  H.  Sampson,  H.  L.  Zhang,  and  C.  J.  Fontes,  At  Data  Nucl.  Data  Tables  48, 25  (1991). 


112 


19.  A.  Dasgupta  and  K.  G.  Whitney,  At  Data  Nucl.  Data  Tables  58, 77  (1994). 

20.  Atomic  stracturc  code  of  R.  D.  Cowan  foUowing  his  book  The  Theory  of  Atomic  Structure  and 
Spectra  (University  of  California  Press,  Berkeley,  CA,  1981). 

21.  H.  L.  Zhang,  D.  H.  Sampson,  R.  H.  Clark,  and  J.  B.  Mann,  At  Data  Nucl.  Data  Tables  37, 17 
(1987). 

22.  P.  Hagelstein  and  R.  K.  Jung,  At  Data  Nucl.  Data  Tables  37, 121  (1987). 

23.  L.  B.  Golden,  D.  H.  Sampson,  and  K.  Omidvar,  J.  Phys.  B  11, 3235  (1978). 

24.  D.  L.  Moores,  L.  B.  Golden,  and  D.  H.  Sampson,  J.  Phys.  B  13, 385  (1980). 

25.  R.  Bartiromo,  F.  Bombarda,  and  R.  Giannella,  Phys.  Rev.  A  32, 531  (1985). 

26.  A.  H.  Gabriel  and  K.  J.  H.  Phillips,  Mon.  Not  R.  Astr.  Soc.  189, 319  (1979). 


113 


TABLE  1.  IDENTIFICATION  OF  NE-LIKE  LINES  AND  RATIOS 


Transitions 

EnerEy(eV) 

Ri 

[2p3/23t^5/2]i  — ^  ^^0 

1574.9 

R2 

[2p3/2  3^3/2]!  — »  ^Sq 

1613.5 

Si 

[2^3/2 3«^5/23t^5/2] 7/2  - ^  [3ti5/2l 

1557.94 

82 

[2ps/23<^S/23t^5/2]5/2  - ^  P^5/2] 

1562.17 

S3 

[2j5l/23tis/23d5/2]5/2  - *■  [^^5/2] 

1598.96 

S4 

l^Pl/2^d$/2^ds/2]7/2  - - ^  [3^^572] 

Ss 

[2px/23cf3/23d3/2l3/2  [3t^3/2] 

1602.39 

and  S2  are  satellites  to  Ri 
83,84  and  8s  are  satellites  to  R2 

£  will  use  the  designatioiis  A  and  B  for  ratios  of  power  emitted  by  Si  and  82  to 
[  respectively  and  C,  D  and  E  for  the  same  ratios  for  lines  S3,  84  and  Ss  to  R2 
respectively. 
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t2  3S  STATES 
36  3P  STATES 
60  3D  STATES 


30  3S  STATES 
90  3P  STATES 
150  3D  STATES 


Fig.  13 
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Vin.  CALCULATING  MEASURABLE  X-RAY  QUANTITIES  FROM  SIMULATED  PRS 


In  recent  ycare,  significant  progress  has  been  made  towards  developing  the  tools  necessary 
to  determine  the  plasma  conditions  generated  in  imploding  z-pinches.  Progress  has  also  been 
made  in  improving  MHD  simulations  of  PRS  implosions,  which  improve  the  effectiveness  of  the 
experiments  by  providing  an  inexpensive  way  to  explore  the  consequences  of  varying  experimental 
parameters.  For  the  most  productive  interacdon,  identical  experimental  and  simulation  quantities 
must  be  compared.  To  do  this,  the  available  experimental  measurements  must  be  analyzed  to  infer 
plasma  conditions,  and  the  detailed  theoretical  calculations  must  be  meaningfuUy  inteipreted  usmg 
the  same  methods  and  procedures  that  were  used  to  analyze  the  experiments.  This  section  descnbes 
a  procedure  for  processing  MHD  calculation  results  to  produce  Bteied  pinhole  photographs  and 

spadaUy  resolved  spectral  data,  as  weU  as  time-resolved  power  outputs  in  analogy  to  experimentally 

observed  quantities.  These  quantities  can  then  be  readUy  and  meaningfuUy  compared  with 

experimental  measurements. 


Analysts  off  Experimental  Data 

Experimental  measurements  of  PRS  implosions  yield  a  wealth  of  data,  including  voltage  and 
current  traces,  pinhole  photographs,  yield  measurements  from  calorimeters  and  diodes  (which  may 
be  filtered  to  restrict  the  frequency  band),  and  quantitative  spectra  (energy  emission  m  spectral 
lines).  In  analyzing  this  data,  the  challenge  is  to  infer  from  these  instrumental  readings  the  values 
of  the  physical  quantities  which  are  most  likely  to  have  produced  the  observed  yield,  spectrum  and 

pulsewidth  of  the  implosion. 

Hie  NRL  Radiation  Hydrodynamics  Branch  has  developed  an  analysis  proce  ure  at  uses 
calculations  of  average  x-ray  power  and  plasma  diameter,  together  with  x-ray  spectrosc^ic 
measurements,  to  infer  the  average  ion  density  and  electron  temperature  m  the  nnploston  This 
procedure  has  been  successfuUy  used  to  optimize  aluminum  wire-array  implosions.  It  butids  on 
a  comprehensive  theoretical  atomic  model,  but  also  requires  a  careful  and  systematic  treatment 
of  the  experimental  data>  in  order  to  determine  useful  values  for  bulk  parameters  hke  average 
power  and  diameter.  Ihe  procedure  derives  time  and  space  averaged  results  for  the  radiatmg 
region,  and  has  been  used  to  diagnose  time-averaged  measurements.*  It  is  important  to  note  that 
it  can  be  used  with  time-resolved  measurements  as  well,  as  long  as  the  experimental  timescale 
is  long  enough  to  establish  coUisional-radiative  equilibrium  (CRE);  this  condition  is  typtcally 
wen  satisfied  on  the  nanosecond  scales  typical  of  present-day  experiments  and  of  htgh-qnahty 
time-resolved  measurements. 
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Analysis  of  Simulation  Data 

MHD  simulations  involve  detailed  calculation  of  plasma  parameters  at  each  point  in  space 
and  time.  The  plasma  parameters  include  electron  and  ion  density  and  temperature,  ionization 
state,  and  radiated  energy  at  all  important  energies  (line  and  continuum).  This  large  amount  of 
information  must  then  be  reduced  to  that  which  is  observed  in  the  experiment  in  order  to  make 
meaningful  comparisons  with  experiment  The  reduction  process  will  necessarily  depend  on  the 
nature  of  the  observation;  for  example,  the  plasma  diameter  will  be  different  for  K-shell  photons 
than  for  L-shell  photons.  For  this  reduction,  we  are  primarily  interested  in  obtaining  the  quantities 
used  in  the  experimental  data  analysis,  namely  the  radiated  x-ray  power  and  spectrum  (at  least 
including  the  Lyman  and  Helium  series),  and  the  plasma  diameter. 

The  main  objective  here  is  to  conduct  the  same  analysis  with  simulation  data  as  was  done  with 
the  experimental  data,  and  compare  the  two  as  a  function  of  time.  Accordingly,  we  used  a  set  of 
aluminum  wire-array  implosions  on  Saturn  with  a  fairly  high  x-ray  yield  and  a  very  symmetric 
implosion  profile,  and  modeled  these  with  a  1-D  MHD  simulation.  The  MHD  code®  has  been 
used  extensively  in  the  past,  and  includes  a  comprehensive  atomic  radiation  model.  It  regards 
the  cylindrically  symmetric  plasma  as  being  divided  into  30  zones,  which  are  annular  shells.  The 
MHD  run  used  the  same  circuit  parameters  as  the  experimental  shot,  and  obtained  a  very  similar 
pulse  profile  and  x-ray  power  (cf.  Fig.  (1)  and  Fig.  (2)).  (The  most  notable  difference  is  the  absence 
of  a  “bounce”  after  stagnation,  due  in  part  to  the  arbitrary  zeroing  of  current  after  stagnation). 
Other  data  measured  or  inferred  from  experiment  are  the  calculated  pinch  radius  r(f),  and  the 
aluminum  K-shell  spectra.  The  experimental  values  for  the  time- varying  electron  temperature  and 
ion  density,  as  determined  by  the  NRL  procedure  discussed  above,  are  shown  in  Fig.  (3).  The 
range  of  MHD-determined  temperature  and  density  can  be  seen  in  Fig.  (4),  which  show  these  as  a 
function  of  time  in  each  zone. 

It  is  evident  from  these  figures  that  the  complete  set  of  simulation  temperatures  and  densities 
can  vary  widely  from  the  experimental  values.  The  diagnostic  significance  of  this  variation 
depends  largely  on  the  radiative  importance  of  the  region;  zones  with  very  low  mass  or  which  are 
obscured  from  the  detector  by  intervening  optically  thick  plasma  might  have  a  significantly  different 
temperature  or  density  from  the  rest  of  the  plasma,  without  affecting  the  radiation  measurements. 
These  variations  might  nevertheless  significantly  change  the  implosion  dynamics,  and  a  future 
application  of  the  present  work  will  be  to  find  ways  to  detect  such  otherwise  hidden  plasma 

variations. 

Theoretical  Generation  of  Pinhole  Kiotographs 

A  pinhole  photograph  is  a  direct  representation  of  the  radiative  profile  of  the  source.  The 
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pinhole  image  at  a  point  is  the  sum  of  the  radiative  emission  and  absorption  in  the  plasma  along 
the  ray  that  connects  the  image  point  and  the  source.  This  is  depicted  in  Fig.  (5),  which  also  gives 
some  of  the  geometrical  relations  characterizing  the  problem  (discussed  below).  The  darkness  of 
a  point  in  the  pinhole  image  should  be  proportional  to  the  energy  that  reaches  that  pomL  For  a 
time-resolved  (or  a  streak)  image,  this  can  be  taken  as  proportional  to  the  instantaneous  radiative 
power  at  the  time  of  the  image,  but  for  a  time-integrated  image,  the  power  must  be  mtegrated  to 

obtain  the  energy. 

ITie  intensity  along  each  ray  in  Fig.  (5)  is  given  by  the  equation  of  transfer. 


_  T  C 

-7 -  —  Iu~ 

dT„ 


(1) 


where  h  is  the  specific  intensity  at  frequency (units,  e.g.,  W  cm-),  and  S.  is  the  source  function, 
T^,  is  the  optical  depth  along  the  ray: 

T>a(a)=  f  i^{a)ds  (2) 

J  So 

and  is  the  absorption  coefficient  (units,  e.g.,  cm-').  The  sign  in  Eq.  (1)  is  chosen  such  that 
the  optical  depth  increases  from  the  observer  towards  the  radiation  source.  We  wish  to  solve  this 
equation  for  the  specific  intensity  at  the  observer.  The  formal  solution  of  the  equation  of  transfer 

at  the  near  point  ro  is  ^ 

=  f  S(r')e-'dr'  +  /(r)e'-"  (3) 

Jto 

Die  boundary  condition  is  that  the  specific  intensity  in  the  direction  of  the  observer  vanishes  at 
the  far  edge  of  the  plasma.  If  the  outermost  zone’s  optical  depth  along  the  ray  is  At,  the  specific 

intensity  at  the  inner  edge  of  this  zone  is  then 

Joue  =  [1  -  .  W 

The  intensity  at  each  successive  inner  zone  along  the  ray  is  readily  found  by  extending  this 
procedure,  using  in  Eq.  (3)  to  get  the  specific  intensity  at  the  next  inner  zone,  and  so  on;  the 
specific  intensity  at  the  plasma  edge  facing  the  observer  gives  the  radiative  energy  produemg  the 


pinhole  image. 

Hie  source  function  S.  and  absorption  coefficient  are  detennined  by  the  plasma  conditions 
in  the  implosion,  and  these  functions  are  in  fact  used  for  energy  balance  in  the  MHD  model.  The 
MHD  radiation  ffansport,  however,  is  concerned  primarily  wi±  eneigy  balance  (radiation  to  fuU  4ir 
steradians);  in  this  case,  the  probability-ofrescape  method'  is  effective,  and  is  extremely  efficient. 
For  pinhole  photographs,  on  the  other  hand,  the  main  concern  is  ray  tracing  and  energy  transport 
to  the  recorder,  and  the  procedure  outlined  above  must  be  used,  with  multifrequency  transport  of 
each  spectral  line  along  each  ray. 


TTie  frequency-dependent  intensity  on  the  edge  of  the  plasma  can  be  used  to  construct  synthetic 
pinhole  images  of  the  evolving  implosion  in  any  wavelength  band.  These  can  be  used  to  detenmne 

the  plasma  diameter  and  spectral  properties  at  the  observer. 

For  experimental  observations,  filters  are  commonly  used  to  restrict  the  spectrum  of  the 
radiation  generating  the  pinhole  photograph.  The  filter  attenuates  the  incident  radiation  along 
each  ray  in  proportion  to  the  filter  thickness  ao  and  the  absorption  function  ni  charactenstic  of  the 


FUtcrs  can  be  used  singly  or  in  combination  to  higUight  features  of  interest  Use  of  experimental 
filter  absorption  functions  is  important  in  processing  theoretical  results  because  the  filters  can 
significantly  alter  the  incident  spectrum,  affecting  the  observed  size  mid  motion  of  the  implodmg 
plasma  as  weU  as  possibly  changing  spectral  features  such  as  which  lines  are  most  promment  In 
the  present  calculations,  however,  the  filter  function  has  been  set  to  unity,  in  order  to  focus  on 
opacity  effects  in  the  imploding  plasma. 

An  example  synthetic  pinhole  image  is  shown  in  Fig.  (6)  for  the  Saturn  aluminum  wire  array 
implosions  chosen  for  comparison.  Tlie  effective  pinch  diameter  is  plotted  as  a  function  of  time  in 
Fig.  (7),  where  is  can  be  seen  that  it  is  at  greatest  compression  about  a  factor  of  2  larger  than  the 
experimentaUy-observed  pinch  diameter,  also  shown.  For  this  presentation,  a  simplified  version 
of  the  pinhole  photograph  procedure  was  used,  which  assumed  that  the  radiative  emission  from 
each  zone  was  isotropic,  and  that  the  net  emission  at  any  radial  position  is  given  by  the  total  power 
absorbed  or  emitted  in  that  zone.  This  simplified  model  does  not  include  the  geometrical  factors 
that  are  found  in  the  detailed  ray-tracing  procedure  described  above,  but  the  predicted  diameter 
is  nonetheless  close  to  that  inferred  from  experiment  (and  to  the  mass-averaged  diameter).  These 
simpUfied  results  are  intended  as  an  computationally  inexpensive  illustration  of  the  method;  later 
calculations  will  use  the  more  accurate  procedure. 


Table  L  Preliminary  results  from  analysis  of  simulation  data 


Time 

a  ratio 

Radius 

Power 

T. 

Ui 

-1.2  ns 

0.85 

2mm 

3.2  TW 

350  eV 

2.2x10*°  cm“* 

0.3  ns 

1.40 

2mm 

3.8  TW 

550  eV 

1.8xl0*°cm“® 

2.4  ns 

131 

3mm 

2.8  TW 

550  eV 

9  X 10^°  cm“* 

4.8  ns 

0.71 

4mm 

1.4  TW 

400  eV 

7  X 10^°  cm"° 

If  the  simulation  results  are  analyzed  using  the  same  procedure^  as  the  experimental 
measurements,  with  the  pinch  diameter  given  as  above,  x-ray  power  obtained  directly  from  the 
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simulation,  and  spectral  line  enetpes  (Ly  «  and  He  n)  given  by  the  net  pmch  “2”': 

™lues  can  be  found  for  the  effective  electron  temperature  and  ion  density  at  times  nem  the  pe 
the  radiation.  These  results  are  shown  in  Table  1,  and  are  generally  comparable  to  those  inferred 

from  experimeiit  (shown  in  Fig.  (3)). 


We  have  briefly  described  a  procedure  for  expressing  theoretical  data  in  the  same  vocabulary 
as  experimental  observation.  Future  work  will  capitalize  on  the  theoretical  capabdity  to  construct 

“what-if  ’  scenarios.  Goals: 

1,  Use  filter  and  device  response  functions  in  the  calculation. 

2  Identify  what  radial  parts  of  the  plasma  dominate  in  producing  each  spectral  fine.  Frequency 

diaenosis  of  the  radial  (temperature,  density)  profile 
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Figure  1.  Experimentally  observed  x-ray  power  from  a  single  implosion  on  Saturn  (shot  2095),  as 
measured  by  a  PCD.  The  experiment  used  an  array  90  aluminum  wires  with  total  mass  300  pg/cm. 
Total  reported  x-ray  yield  was  53  kJ. 
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Figure  2.  TheoreticaUy  (MHD)  predicted  x-ray  power  for  a  wire-^y  implosion 
with  the  same  circuit  and  load  parameters  as  were  used  in  the  implosion  or  Fig.  (1) 
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Figure  3.  Experimentally-determined  electron  temperature  and  ion  density  for  Saturn  alummum  w^e 
array  implostons  (90  wires),  obtained  according  to  the  procedure  in  Ref.3  The  pomts  represent  data 
for  three  different  shots  (2085, 2094  and  2095),  aU  of  which  used  identical  parameters  and  generate 
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Basic  Ray  tracing  geometry 
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Fisure  7  Pinch  diameter,  as  measured  experimentally  (top)  and  as  determined  from  MITO- 

^tive  output  (bottom).  Experimental  diamter  ^  dete^ed  by  ftog  a  Gaussian 
to  the  axially  averaged  profile;  symbols  correspond  to  the  3  shots  modeled,  as  m  F  g. 


ta  last  year’s  final  report,  a  series  of  experiments  carried  out  at  Physics  International’s  Double 
eagle  facility  using  Mg-coated  A1  wires  was  described.  Initial  analysis  of  the  data  available 
indicated  that  the  K  shell  radiation  was  enhanced  by  the  use  of  this  near-Z  element  mixture.  Since 
that  Ume,  more  data  has  become  available  and  the  analysis  was  finalized.  None  of  the  preliminary 
conclusions  have  changed.  n,e  K  shell  yield  increased  from  28-50%  for  an  80:20  Al:Mg  mixmre 
compared  to  pure  A1  loads.  Furthermore,  the  Mg  was  found  to  be  hotter  than  the  Al,  600  eV,  as 
opposed  to  Al’s  300  eV.  Time  resolved  data  supports  the  foUowing  picture:  the  Mg  coat  boils  off 

-  i. — ™ ~ 

behind,  and  evenmally  arrives  near  axis  where  it  becomes  half  as  hot  as  the  Mg  that  preceded  tt. 


Webelieve  that  these  results  reflect  both  opacity  and  temperatureeffects.  The  theory  of  opacity 

eritancement  of  mixed  loads  was  presented  a  decade  ago  in  Ref.  1.  When  some  of  the  Al  wire  is, 

in  effect,  shaved  off  and  replaced  with  Mg,  the  two  dominant  He-  and  H-like  1-2  lines  of  K  shell  Al 
are  somewhat  less  optically  thick,  20%  less  in  the  most  interesting  case  of  the  80:20  Al:Mg  mixture. 
Furthermore,  the  strong  analogous  lines  of  Mg  are  subject  to  much  less  opacity  attenuation  than 
those  of  Al.  since  Mg  represents  only  one-fifth  of  the  load.  The  result  is  an  increase  of  the  K  shell 
power  even  if  the  temperature  of  the  elemental  constituents  is  the  same.  The  extra  heating  of  the 
Mg  enhanced  the  effect.  This  data,  representing  an  important  advance  in  out  understandmg  of  Z 
pinches,  has  now  been  published  in  the  Physical  Review  (Ref.  2). 


This  very  encouiging  experimental  result  with  two  specific  elements  represents  the 
confinnation  of  an  effect  predicted  a  decade  ago  in  Ref.  1.  However,  there  is  no  n  priori 
reason  to  believe  that  Al  and  Mg  are  necessarily  the  optimal  pair  of  elements  or  that  more  than  two 
elements  might  not  do  an  even  better  job  of  yield  enhancement.  Note  that  there  is  another  DNA 
goal  which  this  effect  serves  well:  fidelity  enhancement.  The  use  of  mixtures  provides  x-rays  m 
different  regions  of  the  spectrum  rather  than  concentrating  them  in  one  narrow  region  of  K  sheU  Ime 

The  fact  that  there  are  SO  many  easily  adjustable  parameter,:  the  number  of  elements 

used,  their  atomic  numbers,  and  their  proportions  in  the  load  gives  reason  for  optimism  that  very 


useful  combinations  will  be  found  and  become  a  staple  ofDNA’s  simulation  capability. 

To  get  closer  to  a  determination  of  what  would  be  optunal,  a  series  of  calculations  was  done 
using  all  of  the  elements  of  atomic  numbers  between  1 1  and  14:  Na,  Mg,  Al,  and  Si,  respectively. 
Temperatures  of  100  to  1000  eV  were  examined,  covering  most  of  the  range  in  which  these  elements 
are  stripped  to  the  K  shell.  Fig.  1  shows  the  fraction  of  each  element  which  exists  in  the  H-like  or 
He-like  stages  as  a  function  of  temperature  for  typical  pinch  conditions.  This  calculation  displays 
quantitatively  the  need  for  increasing  temperature  with  increasing  atomic  number.  The  K  shell 
ionization  fraction  for  each  element  is  a  major  determinant  of  its  ability  to  radiate  copious  K  shell 

emission. 

Two  other  factors  are  also  important  in  determining  the  emitted  K  shell  power  for  a  given 
element.  First  is  the  collisional  excitation  rate  for  the  levels  which  radiate  the  lines.  This  rate 
generally  increases  with  temperature.  The  other  important  physical  process  is  self-absorption  of 
the  line  radiation  due  to  the  generally  high  opacity  of  the  strong  lines.  The  result  of  these  latter  two 
factors  is  that  the  radiated  power  tends  to  increase  even  when  ionization  reduces  the  fraction  of  the 
radiating  ion.  As  ionization  proceeds,  the  optical  depth  of  the  strong  lines  is  reduced  since  there  are 
fewer  absorbers  in  the  line  of  sight.  The  increased  collisional  creation  of  photons  further  augments 
the  emission.  The  result  for  these  four  elements  is  displayed  in  Fig.  2,  which  presents  the  K  shell 
radiated  power  as  a  function  of  temperature  up  to  a  keV.  In  each  case,  a  rapid  initial  rise  with 
temperature  is  followed  by  a  gentler  increase  at  higher  temperatures.  Note  that,  for  temperatures 
in  excess  of  400  eV,  the  higher  the  atomic  number,  the  greater  the  radiated  power.  A  key  effect 
contributing  to  this  result  is  the  atomic  number  (Z)  scaling  of  the  optical  depth  of  a  given  transition. 
Such  scaling  depends  upon  the  line  Doppler  width  and  some  fundamental  atomic  parameters,  and 
shows  that  the  optical  depth  of  a  specific  line  at  a  given  density  varies  as  Z~^ ,  i.e.,  a  sharp  reduction 
in  opacity  with  increasing  atomic  number.  Therefore,  if  Si  can  be  stripped  to  the  K  shell,  it  is  a 
better  radiator  than,  for  instance.  Mg.  For  the  pinch  conditions  found  in  the  Al:Mg  experiments 
performed  at  PI,  we  find  that  had  the  Mg  coat  been  Si,  the  K  sheU  yield  would  have  doubled  rather 
than  increasing  by  up  to  50%.  A  lower  atomic  number  element  such  as  Na  would  be  useful  if  the 
load  can  be  designed  to  place  it  in  the  cooler  outer  region  of  the  pinch.  It  is  clear  that  both  yield 
and  fidelity  can  be  improved  by  further  exploration  of  these  effects  in  well-considered  experiments. 
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One  of  the  more  editing  aspects  of  the  mixture  concept  is  that  higher  atomic  number  elements 
hardct  photons  ate  even  more  susceptible  to  positive  results  than  the  elements  already 

successfully  demonstrated. 
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FIG.  1.  Sums  of  the  H-  and  He-like  fractions  are  plotted  as  a  function  of  electron  temperature  for 
the  indicated  plasma  conditions  for  Na,  Mg,  Al,  and  Si. 
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FIG.  2.  K  shell  line  power  is  shown  vs.  electron  temperature  for  Na,  Mg,  Al,  and  Si,  for  the 
indicated  pinch  conditions. 


